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Abstract 

This article is devoted to the well-posedness of the stochastic compressible Navier Stokes equations. We 
establish the global existence of an appropriate class of weak solutions emanating from large inital data, set 
within a bounded domain. The stochastic forcing is of multiplicative type, white in time and colored in space. 
Energy methods are used to merge techniques of P.L. Lions for the deterministic, compressible system with 
the theory of martingale solutions to the incompressible, stochastic system. Namely, we develop stochastic 
analogues of the weak compactness program of Lions, and use them to implement a martingale method. 
The existence proof involves four layers of approximating schemes. We combine the three layer scheme of 
Feiresil/Novotny/Petzeltova for the deterministic, compressible system with a time splitting method used 
by Berthelin/Vovellc for the one dimensional stochastic compressible Euler equations. 


1 Introduction 

This paper is devoted to the analysis of the initial boundary value problem for the compressible isentropic 
stochastic Navier Stokes equations. This system of stochastic partial differential equations governs the evolution 
of a viscous, compressible fluid(or gas) subject to random perturbations by noise. The macroscopic state of the 
fluid is described by a pair (p, it) consisting of the scalar, nonnegative density p and an M. d valued velocity field 
u. In the isentropic case, the system is written as 

{ d t p + div (pu) = 0 

dt(pu) + div (pu ® u) + VP(p) — 2/iA u — AV div u = pok(p, pu , x)$k (1.1) 

0(0), 00(0)) = Oo, mo) 

where p, A > 0 are positive viscosity coefficients and P(p) is the macroscopic pressure. The collection {/3k}%L 1 
consists of independent, R valued Brownian motions and {<Jk}kLi is a collection of R d valued noise coefficients. 
We fix a large time T and pose the system ED on (0, T] x D , supplemented with a Dirichlet boundary condition 
for the velocity. 

The main contribution of this article is as follows: we introduce a notion of global weak solutions to ED 
and provide an existence proof under suitable hypotheses on the pressure, the data, and the noise. 

1.1 Hypotheses and Relevant Literature 

Let us begin with our assumptions about the pressure law, the initial data, and the noise coefficients. As is 
common in the literature for the deterministic equation, we impose an equation of state for the pressure law. 

Hypothesis 1.1. The pressure law takes the form P(p) = p 1 ■ In dimensions d = 2,3 we impose 7 > f,f 
respectively, and for d > 3 we assume 7 > |. 

The initial data are a finite energy density/momentum pair satisfying a compatibiility criterion. 

Hypothesis 1.2. The initial data (po,rno) are deterministic. The initial density po is nonnegative and the 
following compatibility condition holds 

m 0 l {po=O }=0. (1.2) 


1 


Moreover, the initial energy is finite: 
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dx < oo. 


(1.3) 


We make two types of assumptions about the noise coefficients. The first is a continuity hypothesis for each 
k fixed and the second is a trace class type summability condition for the collection {crfe}^ =1 . 

Hypothesis 1.3. For each k, the coefficient a k : R + x x D —> is bounded and continuous. Moreover, 
the following uniform lipschitz condition holds: there exists a constant C k such that for each mi, m 2 £ 


sup \a k (p,mi,x) - a k (p,m 2 ,x)\ < C k \mi - m 2 \. 

(p,a?)£M+ xD 


(1.4) 


In addition, for each pi, P 2 £ 


sup \a k (pi,m,x) - a k {p 2 ,m,x)\ < C k \pi - p 2 \. 

(m,x)GK d X D 


(1.5) 


Hypothesis 1.4. The sequence of coefficients {o' k }^L 1 satisfy the summability relation 
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where we denote for p > 1 

Wk\= \snp\a k (p,m,-)\\ P . 

P ’ P ^ 

A few remarks are in order. A simple, but important point is that the initial data are described by a density/ 
momentum pair (po,mo) rather than a density/velocity pair. As a result, the initial velocity is undefined in the 
vaccum regions. Next we note that the assumptions on 7 imposed in Hypothesis 1 1.1 1 are slightly worse than in 
the original work of Lions m in dimension 3, in that we ask for a strict inequality 7 > |. This simplifies the 
proof but is not necessary. Using the tools developed by Feiresil in 0, m one can weaken the constraints to 
simply 7 > f in all dimensions. 

Regarding Hypothesis 11.41 a simple example to have in mind is a noise of the form pW where W [x, t) = 
X)fc°=i V^k e k{x) for some orthonormal basis {efc }^ =1 of [Lff\ d where 

OO 

YMe k \ 2 -/ < 00 . 

“ L ■y ~ 1 

k=l 

In this case, a pathwise approach to the solvability of system (O) is possible, as shown in m- 

The starting point for our analysis is a formal energy identity (see Section ^. 31) . In the stochastic, compressible 
framework, the kinetic and potential energy are random processes which fluctutate due to noise and grow 
according to an Ito correction term. To close an estimate on their moments(in w) and obtain an a priori bound 
for the SPDE, one is lead to the trace class type summability condition for the coefficients {u k } k L 1 (Hypothesis 
11.41) . This leads one to suspect that for reasonably smooth noise coefficients satisfying Hypothesis 11.41 finite 
energy weak solutions to cni) should exist. 

Unfortunately, to construct these solutions, a simple Galerkin approximation procedure is not enough. 
Instead, one requires several layers of approximating stochastic PDE’s. Within each layer, one establishes 
a priori bounds and implements a compactness method. Broadly speaking there are two principle steps, an 
application of a generalized version of the Skorohod theorem (Theorem 1A. 121) to obtain a form of compactness 
and an alternative to the martingale representation theorem (Lemma 1A. 151) to pass limits in the stochastic 
integrals. However, there are a number of subtleties in the implementation due to the limited control that the 
a priori bounds provide on the density. 

The literature devoted to the the deterministic, compressible system is extensive; the most fundamental 
for our work are the results of Lions [21j and Feiresil/Novotny/Petzeltova [12] . These works provide the main 
inspiration for the subtle analysis of the density and the particular approximating schemes, respectively. 
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There has also been an intensive study of the incompressible, stochastic Navier Stokes equations. The 
works most relevant to our paper concern the construction of weak, martingale solutions. See for instance 
Mm m- There are also a few articles concerning the non-homogenous, incompressible system; see |23] 
and [7]. However, the literature concerning the stochastic, compressible system is rather scarce. Some results 
are available in dimension one, see [24]. The most relevant for our analysis is the work of Berthelin/Vovelle [3] 
on the one dimensional compressible Euler equations. This paper inspired the time splitting scheme used in the 
lowest level of our approximations. 

The paper nn studies the compressible Navier Stokes equations driven by a forcing of the form pW. In this 
special case, one can change variables and work pathwise in to. This technique is not generally available for 
other types of multiplicative noise. 

The only existing literature on the compressible, stochastic Navier Stokes equations that is comparable to 
this work is a very interesting preprint of Breit/Hofmanova [3] . We became aware of these results during the late 
stages of the write up of this work. We emphasize that our results were independently concieved and obtained. 
There are a number of similarities between their work and ours, but our hypotheses on the noise do not overlap, 
so neither result implies the other. As a result, there are also several differences in the approximating schemes. 
Finally, we should remark that our work is set on a bounded domain, rather than the torus. 

1.2 Notion of Weak Solution and Statement of the Main Result 

In this article, we refer to (fl, T, IP, {.Ftl^g, as a stochastic basis provided (f2, J - , P) is a probability 

space endowed with a collection {fJ k Yk—i of one dimensional, independent Brownian motions adapted to {fFt}t=o 
(see Appendix A for more discussion). 

The solutions we construct are weak in both the analytic and the probabilistic sense. Namely, a solution to 
m is a pair (p, u ) satisfying the continuity and momentum equations in the analytically weak sense, relative 
to a stochastic basis of our choice. More precisely, weak solutions are defined as follows: 

Definition 1.5. A pair ( p,u ) is a weak solution to the stochastic compressible Navier Stokes equations (11.11) 
provided there exists a stochastic basis (O, Jq {.F 4 }^, P, {/3fc} fc>1 ) such that 

2~l 

1. The pair (p, pu ) : LI x [0, T] —> [Lf x Lf +1 is an {F^j^-g progressively measurable stochastic process with 
P a.s. continuous sample paths. The velocity u £ L 2 (LI x [0,T];ffg. I ,) is an equivalence class o/jF^l^lg 
progressively measurable Hf valued processes. 


2. Continuity Equation: For all <f> £ Cf°(D) and t £ [0, T], the following equality holds P a.s. 

/ p{t)4>dx = / po<j)dx + / / pu ■ V <f>dxds. 

J D J D Jo J D 

3. Momentum Equation: For all </> £ [C£°(D)] d and t £ [0,T], the following equality holds P a.s. 

/ pu(t) ■ <f>dx = / mo ■(/)+ / [pu ®u— 2/iVrt — A div ul] : \7(f> + p 7 div (f> dxds 

Jd Jd Jo Jd 


°° ft r- 

+ V / / pcr k (p, pu,x) ■ 4>dxd/3 k (s). 
k=1 J o Jd 


4- The following energy estimate holds: for all p > 1 
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We can now give a precise statement of our main result: 

Theorem 1.6. Let (po, m o) be an initial density/momentum pair satisfying Hvvothesis \1.2\ Suppose the pressure 
P(p ) satisfies Hypothesis \1.1\ Let {a k }k >l be a collection of noise coefficients satisfying Hypotheses 11.11 and 

f T3L 

Then there exists a weak solution (p, u) to the stochastic compressible Navier Stokes system (11.11) in the 
sense of Definition \1.5l 
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1.3 Outline of the Proof and Main Difficulties 

Let us proceed to an overview of the proof of Theorem 11.61 The construction of weak solutions involves four 
layers of approximating schemes; labelled (from the lowest to the highest level) r, n , e, and 8. 

In Section [31 we prove existence for the lowest level of our approximating scheme. Definition 13.31 introduces 
the notion of a r layer approximation to the compressible, stochastic Navier Stokes system m- The main result 
of the section, Theorem 13.11 shows that for each fixed (n, e, 8), one can construct a sequence {(/5 r , w, T )} T >o of r 
layer approximations which are controlled uniformly (in all parameters, though we will only indicate dependence 
on one parameter at a time) in the sense of the energy estimate (13.21) . 

The proof is based on a time splitting method used by Berthelin/Vovelle in [3]. Half of the time, the 
determinstic system evolves, and the stochastic forcing is neglected. The other half of the time, the density 
is frozen, and the system evolves only through the noise. The evolution is sped up appropriately so there is 
consistency when the time splitting parameter is sent to zero (in Section [¥]). The main tools of this section are 
Propositions 13.61 and 15771 which use classical fixed point arguments to obtain a basic existence result for both 
the deterministic and the stochastic systems, respectively. The solutions are sufficiently regular so that Ito’s 
formula may be applied, and the formal estimates of Section [2731 can be justified. This leads to uniform bounds. 

Section[4]is devoted to proving Theorem l4.ll an existence result for the second layer. Definition l3.3l introduces 
the notion of an n layer approximation CH), and the goal is to construct a sequence {(pn,w n )}^ =1 of these 
which obey the uniform bounds (14.11) . 

The proof uses the existence theory at the r layer together with a compactness method. For each fixed 
n > 1, we apply Theorem 13.11 and find a sequence {(/5 Tjn , m Tj „)} t >o of r layer approximations obeying the 
uniform bounds (13.21) . Our first step is to establish Proposition [4731 which yields a limit point ( p ni u r ), together 
with a new sequence {(p Ti „, m Ti „)} t >o of r layer approximations with improved compactness properties. A 
caveat is that the new sequence is defined relative to a new probability space. However, there exists a sequence 
{T t } t> o of measure preserving transformations (referred to by the author as recovery maps), which link the 
new and the old sequence by composition. These mappings allow us to preserve information as we change 
probability spaces; in particular, ensuring that the new sequence solves the same equations, obeys the same 
uniform bounds and is generally unaltered in any of its arguments besides oj. Our main tool is Theorem IA.121 
a generalization of the classical result of Skorohod. It merges two recent extensions of the theorem; the first, 
due to Jakubowksi [16] . permits random variables on a class of topological spaces and the other, due to by 
Vaart/Wellner [26], provides the recovery maps. 

A distinctive feature of the r layer is that one requires stronger convergence in time than in other layers of 
the scheme. This is required to compensate for the high frequency switching between the two types of evolution 
when t is small. This leads to some subtleties in the proof of tightness Lemma 14.41 To obtain the necessary 
Holder estimate in time, one needs a probabilistic bootstrapping procedure to deal with the coupling between 
the density and the velocity. 

After Proposition 14.31 is established, it is easy to pass to the limit in the parabolic equation on the new 
probability space ('Lemma 14.51) . To pass to the limit in the momentum equation (Lemma 14.61) . we use a 
martingale method based on an appendix result lA. 151 which provides a convienient characterization of a series of 
one dimensional stochastic integrals. This was developed in [5] as an alternative to the martingale representation 
theorem. This method is used systematically throughout the paper when passing to the limit in the momentum 
equation at each layer. 

Section [5] is devoted to existence Theorem 15.11 the e layer. For each e > 0 fixed, we can apply Theorem 14.II 
to obtain a sequence {(p n ,e, w n ,e)}^i of n layer approximations satisfying the uniform bounds (14.11) . Again, the 
section splits into a compactness step, Proposition 15.31 and an identification step, Lemmas 15.51 and 15.61 In this 
section, the spaces where the tightness Lemma |5.41 are proved become a bit more sophisticated. In particular, 
we must use certain Banach spaces endowed with their weak or weak-* topology. At this stage, and at all later 
compactness steps, the Jakubowski extension of the Skorohod theorem is essential. 

As n —> oo, the most challenging term is e'\/u n \/p n in the momentum equation, which corrects for the 
vanishing viscosity regularization in the energy balance. To treat this difficulty, we adapt a technique of Feireisl 
in Lemma 15.51 and upgrade the convergence of the density. This allows us to use a martingale method again in 
Lemma l5.6l and complete our stability analysis at this layer. Finally, we appeal to standard lower semicontinuity 
arguments to uphold our uniform bounds. 

In Section [6] we build our final approximating scheme, a sequence {(pa, Ui)}i>o of 8 layer approximations to 
EID- Again, the proof still splits broadly into two parts; a compactness step and an identification procedure. 
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However, as e —> 0 we encounter several new difficulties related to the pressure V(pJ + $Pe )• The first is that 
basic energy bounds only provide moment estimates on the L\ x norm of the pressure. To obtain tightness 
of the pressure sequence, we must improve these bounds. In Proposition 16.41 we prove stochastic analogues 
of the integrability gains observed by Lions m- Namely, we show that our weak solutions inherit additional 
integrability from the equation itself. 

An even more serious difficulty is passing to the limit in the pressure, which requires strong convergence of 
the density. The compactness step, Proposition 16.31 is setup in a way that anticipates this problem. We start 
with a preliminary identification step, passing to the limit in the continuity equation ('Lemma 16.6D and also the 
momentum equation fLemma 16.7H . modulo a possible change in the pressure law. 

Next we improve upon our preliminary identification step and work towards the strong convergence of the 
density. In Lemma 16.81 we prove a stochastic analogue of the weak continuity of the effective viscous pressure, 
originally discovered by Lions m- Namely, the weak continuity holds after averaging out the contribution of the 
stochastic integral. In Lemma 16.91 this result is used together with techniques from the theory of renormalized 
solutions of the transport equation [ 8 ] to prove the strong convergence of the density and complete the 8 layer 
existence proof. 

At this stage, we have suceeded in constructing a sequence {{ps, m< 5 )}< 5 >o of 8 layer approximations to the 
stochastic Navier Stokes equations which obey the uniform bounds (16.11) . We now proceed to the prove of our 
main result, Theorem m As usual, we begin with a compactness step, Proposition 17.11 The statement of 
the theorem is rather technical, but natural if one anticipates again the difficulties with the pressure term. In 
Proposition 17.161 we prove another integrability gain, which leads to moment bounds of the norm of the 
density. The assumptions on 7 in Hypothesis 11.11 ensure p > 2. This avoids difficulties with renormalizing the 
transport equation. 

The identification procedure faces a new difficulty at this stage, regarding the nonlinear compositions in the 
multiplicative noise. Namely, at each of the earlier stages in the analysis, we used the parameter 8 to regularize 
the density and momentum before composing with the diffusion coefficient <7k- This was crucial for checking 
the two parts in our key identifiation Lemma IA.15I In the 8 layer existence proof, it allowed us to make a 
preliminary passage to the limit in the equation before proceeding to the proof of the strong convergence of the 
density. In this final step, this is simply not possible. Instead, we can only prove that the momentum process 
minus its drift is a martingale iLemma 17.51) . Nonetheless, we show that this is in fact enough information to 
prove again a stochastic analogue of weak continuity of the effective viscous flux. Namely, Lemma 17.61 uses 
the momentum martingale together with a regularization procedure to establish an averaged Ito product rule, 
which is in turn enough to establish the weak continuity, Lemma 17.71 This is used again to prove convergence 
of the density in Lemma l7.8l and strong convergence of the momentum in Lemma 17.91 Finally, we conclude the 
proof of our main result by passing the limit in the momentum equation, once and for all in Lemma 17.51 

The remainder of the paper is organized as follows. Section [2] contains a preliminary discussion of our 
result, while Sections 131161 are devoted to the proof of Theorem 11.61 Section 12.11 reviews the key ideas from the 
literature; Section 12.21 sets the notation; Section 12.81 shows the formal energy estimates; Section 12.41 discusses 
the hypotheses in more detail. Section [3] establishes existence at the r layer; Section [4] establishes existence at 
the n layer; Section [5] builds solutions to the e layer; Section [ 6 ] builds solutions to the 8 Layer; Section [7] passes 
to the limit in the 8 layer to complete the proof of Theorem 11.61 

2 Preliminaries 

2.1 Background/Literature Review 

The section consists of some high level remarks regarding the existence theory for the deterministic, compressible 
system and the incompressible, stochastic system. Both these theories strongly influenced the methodology in 
this paper, so we choose to review some of the main ideas. The reader may wish to skim this section on a first 
reading since, strictly speaking, we do not quote directly from either of the theories. 

2.1.1 Existence Theory for the Deterministic, Compressible System 

The literature on the deterministic system is extensive, and we will not attempt to give a complete discussion of 
the current status of the field. Instead, we focus on the results that provide the guiding principles for our work. 
The seminal work of P.L. Lions |2Tj initiated a large data global existence theory for finite energy weak solutions. 
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Let us give a very rough outline of the construction. The proof splits into two parts; proving that the solution set 
is weakly compact and constructing several layers of approximating schemes. That is, suppose {{p n , u ra )}^= 1 is 
a sequence of weak solutions (or well chosen approximate solutions) which are uniformly bounded in the natural 
energy space. The strategy is to show that if the initial data are stongly convergent, then the corresponding 
solutions must converge to a solution (p, u ) emanating from the limit point of the data. Since the pressure is a 
nonlinear function of the density, the only feasible way to proceed is by proving that the sequence of densities 
converge strongly to p. However, the continuity equation is driven by too rough of a velocity field 
to provide any control on the densities in a positive Sobolev space. Hence, this is a nontrivial task and the 
basic energy bounds alone are not enough. Nonetheless, Lions found a more subtle mechanism in the nonlinear 
structure that gives compactness. 

To motivate the proof, recall the method for obtaining compactness of the density in the Di Perna/Lions 
[5] theory of the transport equation driven by a “rough” velocity field with bounded divergence. One starts 
with a convienient renormalization(meaning just a smooth function to be applied to the density), for instance 
(3(p) = p 1 2 and renormalizes the equation at the level of both the approximation and the level of the limiting 
solution, known a priori to be a renormalized solution of the transport equation in its own right. If strong 
oscillations in the density sequence are present, the operations of composition with a nonlinear function and 
extraction of a weak limit do not generally commute. However, the renormalized form allows one to track the 
evolution of this “commutator” p 2 — p 2 and a Gronwall argument shows that if compression effects are limited, 
strong convergence of initial densities implies the “commutator” vanishes for all later times. Unfortunately, one 
cannot apply this method directly to the compressible Navier Stokes system because the known a priori bounds 
are not enough to rule out the possibility of extreme compression(or expansion). To proceed, Lions made the 
crucial observation that a sort of “monotoncity miracle” occurs for particular pressure laws, and in some sense, 
it suffices that the so called effective viscous pressure P{p) — (2p + A) div u is “slightly well behaved” (in the sense 
of a certain weak continuity property), even if the divergence of the velocity field alone is potentially unbounded. 
The importance of this quantity had already been observed in a simpler context by D. Serre. Moreover, the 
evolution of this quantity is readily available upon taking the divergence, followed by the inverse laplacian on 
both sides of the momentum equation. By studying this quantity before and after a preliminary passage to the 
limit in the momentum equation, one is able to prove a subtle compactness result, known as the weak continuity 
of the effective viscous pressure, which is just barely enough to complete an analysis of a similar “commutator” 
as in the bounded divergence case, and hence conclude the strong convergence of the density. 

The original work of Lions considered power laws with 7 large enough to ensure that the continuity equation 
could be renormalized, see Hypothesis II.II below. Several years later, Feiresil introduced in [9] some additional 
tools which, combined with Lions general strategy of proof, succeeded in weakening the hypothesis on 7 in 
dimensions two and three, to what seems to be the critical level Q 7 > f ■ This is a nontrival task, since 
for low enough values of 7, one dips below the integrabililty required to classically renormalize the continuity 
equation. More importantly for this paper, with co-authors in 112] Feiresil developed a somewhat simplified(but 
still rather long) approximation scheme, based on a Galerkin appromation for the velocity, a vanishing viscosity 
regularization for the continuity equation, and an artifical pressure regularization. 

2.1.2 Existence Theory for the Stochastic, Incompressible System 

There is also a fairly developed literature concerning the stochastic Navier Stokes equations for incompressible 
fluids, which we will not review in much depth. Naturally, much more is known in dimension two, but the 
existence of weak solutions is known in any dimension. In this regard, the primary inspiration for our work is 
Flandoli’s construction of weak martingale solutions in tni; see also j6| . The main point we wish to emphasize 
is that these solutions are weak in both the analytic and the probabilistic sense. Namely, one is allowed to input 
the probability space where the solution is built, along with a convienient choice of Brownian motions and a 
sufficiently large filtration(potentially larger than the information needed to assess the values of the Brownian 
motions alone). These inputs are referred to as the stochastic basis, and after they have been fixed, one asks 
that the momentum equation holds in an analytically weak sense in space and in the Ito sense in time. 

To understand the virtue of flexibility in the choice of a stochastic basis, recall Leray’s construction of 
weak solutions to the deterministic, incompresible Navier Stokes equations. The key point is that uniform 

1 Below this level, one can just barely give a meaning to the flux term in the momentum equation, and Lions method seems to 

break down. 
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bounds in L 2 (Wf’ 2 ) fl C^H^ 1 ) allow one to apply the Aubin/Lions lemma and obtain strong compactness in 
L 2 (L 2 )( and hence weak stability of the flux term), leading to a straightforward (from a modern point of view) 
weak compactness theory. At a superficial level, in the stochastic case, there is an additional variable w, and the 
possibility of “oscillations” in this variable may block the compactness upgrade from the space/time bounds. 
However, if one is content with only accessing the probability law of the solution, then there is a classical 
fix. Namely, if one can show that the sequence of Galerkin approxmations becomes uniformly concentrated(up 
to a set of very small probability) on L 2 (W~ x ’ 2 ) fl C^ff” 1 ), the Skorohod embedding(for random variables on 
complete separable metric spaces) guarantees the existence of a new sequence of random variables (with the 
same distribution) on the unit interval, along with a limit point, for which the usual L 2 x convergence holds 
pointwise. Essentially, under an appropriate change of variables one is able to convert information that only 
holds on average on the initial probability space, to information that holds in every state of the universe of 
a well chosen probability space. One could visualize this in one dimension by noting that given a sequence 
of bumps sliding back and forth across the unit interval on smaller and smaller measure sets, if we rearrange 
the sequence based on the distribution of mass, one converts the typical counterexample to “weak convergence 
implies pointwise convergence” into a pointwise converging sequence, without altering its probability law. 

2.2 Notation 

Recall that the equation is posed in the space/time domain D x [0, T\. Assume that D is a connected, bounded 
open subset of with smooth boundary. The shorthand notation L q (L p ), L q (W k,p ), Wt’ q (L p ) is used to 
denote the spaces L q ([0, T}\ L P (D)) , L q ([0, T]; W k,p (D)) , W k,q ([0, T]; L P (D)) respectively, where each space 
is understood to be endowed with its strong topology. We will often use the same notation to denote scalar 
functions in L q (L p ) and vector valued functions(with d components) in [L q (L p )] d , but the meaning will always 
be clear from the context. To emphasize when one of the spaces above is endowed with its weak topology, 
we write [L q {L p )] w , [L q {W k,p )\ w . Also, the abbreviation C t ([L p ] w ) denotes the topological space of weakly 
continuous functions / : [0,T] —> L P (D). The space is the closure of the smooth compactly supported 
functions, with respect to the W k,p (D) norm. Moreover, we denote Wq 2 as Hq . Given a probability 

space (f2,.F, P) and a Banach space E , let L p (fl;E) be the collection of equivalence classes of T measurable 
mappings X : Q E such that the p th moment of the E norm is finite. Again, we write (B; E) when 
emphasizing that the space is endowed with its weak topology. To define the sigma algebra generated by 
various random variables, we use a restriction operator r t : C([0, T ]; E) —> C([0, t \; E) which realizes a mapping 
/ : [0, T] —> E as a mapping r t f : [0, t] —► E. The same notation is used for the restriction of an equivalence class 
/ £ E 2 ([0,T];E) to rtf £ E 2 ([0,t];E). We denote A = VA -1 , understood to be well defined on compactly 
supported distributions in R d . The symbol B is reserved for the Bogovski operator, see the remarks preceding 
Lemma |B.5l for the definition of the operator, along with its basic properties. Given two dx d matrices A,B, 
A : B denotes a Frobenius matrix product. The notation A < B denotes inequality up to an insignificant 
constant. The notion of insignificance will be clear from the context. 

2.3 Formal Energy Estimates 

In this section, we present a formal derivation of the basic energy equality for the system uni). The kinetic 
and potential energy dissipate in the usual way, but also fluctuate due to noise and grow according to an Ito 
correction. Nonetheless, we will see that Hypothesis 11.41 ensures that the moments(in oj) are controlled by the 
initial energy. 

As there is no Ito term in the continuity equation, we may apply the ordinary product rulc(in time) to 
dt(pu ), the product rule in space to div(pu ® u), and use the equation for dtp to formally rewrite (11.111 as 

j d t p + div (pu) = 0 

y p[dtu + (u ■ V)u] + V(p 7 ) — 2/iAu — AV div u = pat(p , pu, x)pk- 

Dividing the momentum equation by p and noting the remaining noise term, we see \du\ 2 = \ a k\ ■ To 

derive the energy identity, multiply the momentum equation by u and integrate over D. Note that Ito’s formula 
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gives pdtu ■ u = pdt(^~) — \p\du\ 2 . Using the continuity equation again, 


p[d t u+ (u • V)u] • u = pd t {^-) + pu ■ V(^-) - \ a k 


k =1 


oo 

= dt(~P\u\ 2 ) + div(-p\u\ 2 ) - l afe | 2 - 


fc =1 


Multiplying the continuity equation by 1 and noting div(pit)^ 3 j-p 7 1 = di v(^ 3 jP 7 m) — u ■ V(p 7 ) leads 

to the identity u- V(p 7 ) = 3t(^rj/0 7 ) + div(^ 3 j-p 7 M). Combining these observations, integrating the dissipation 
by parts and using the dirichlet boundary condition gives for all t G [0, T] 

J (^p\ u \ 2 (0 ■+-— j-p 7 (i)^ dx + / J 2p\S7u\ 2 + X(divu) 2 dxds 

-jPoj dx + lt. J d pu-a k (p,pu,x)dxd/3 k (s) (2.1) 


1 l^ol 2 1 

2 Po 


oo „t p 

V / p|ofc(p,pu,x)| 2 da;ds. 
fc=l - / 0 Jd 


To control the moments of the LHS, note that the series of stochastic integrals on the RHS is a martingale with 
quadratic variation(at time T) given by X^fc°=i Iq (Id P u ' a kdx) 2 ds. Hence we may use the Burkholder/Davis/Gundy 
inequality followed by Holder and Hypothesis 11.41 to find for p > 1 
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Maximizing over [0,T], then taking the expectation of the p moment on both sides of the energy identity yields 
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Applying Cauchy’s inequality with the exponent pairs (2y, 2, ^zy) and ( 7 , ^zy) followed by Poincare yields 
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I Jpi 
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2.4 A Few Technical Remarks 

27 

In [21], Lions treats(among other possible assumptions) a forcing of the form pf where / € Lj(Lx -1 ). Hence, 
the spatial exponent in Hypothesis 1 1 .41 is a familiar one. However, note that if we were to treat a time dependent 


































<7fc, the summability criterion would require a norm of a k in Lf (LJL™ n ) in order to close the energy estimates 
as in Section E31 

Next we make a few remarks regarding the notion of weak solution, Definition 11.51 Note that we add noise 
in the spirit of Krylov m by working directly with a series of one dimensional stochastic integrals. Let us 
check that Part [l] of Definition 11.51 combined with Hypotheses 11.3111.41 imply that the series in (11.8|) admits a 
well defined continuous martingale version. For each k, define the process f k by the relation 


fk(s) = / p(x, s)a k {p(x, s),pu(x, s),x) ■ (j>(x)dx. 
J D 


In view of Lemma lA. 141 it suffices to check that fk is an {J^}-f =0 progressively measurable process and 


E E ' 

k= 1 


fk(s)ds 


< oo. 


The desired summability follows from Holder and Hypothesis 11.41 since 


E 1 


[j (yj p{x, s)a k (p(x,s),pu(x, s),x) ■ cj)(x)dxj ds] 


(L° 


E P [|p( s )li-r] ds < \a k \ 2 

L 7 - 1 (L“ ) 


Using the continuity Hypothesis 11.31 one may check the measurability of the following map(by regularization) 


(P,m) G 


LZ x L- 


->■ / p(x)cr k (p(x),m(x),x)dx. 
JD 


Hence, f k inherits progressive measurability from the density/momentum pair (p, pit), in view of PartQ] 

Also, the careful reader may wonder whether there is some ambiguity in the notion of the momentum 
equation [3] regarding measure zero sets(due to the fact that we do not work explicitly with infinite dimensional 
stochastic integrals). However, this is not the case, and one can check that it is equivalent to ask for a universal 
set of full P measure where the weak form holds for all <f> G [Cf°(D)] d simultaneously. This follows from a 
density argument provided one chooses a suitable modification of the stochastic integrals. 

We should also mention that it is unclear(to the author) whether the velocity field u inherits even a weak 
form of continuity in time from the equation itself. This is the reason why we do not ask that the velocity 
u is a stochastic process in the usual sense(unlike the density/momentum pair), and is only identified up to 
equivalence relations in a class of non-anticipating processes. Also note that the dirichlet boundary condition 
is understood in the weak sense. 


3 r Layer Existence 

In this section, we build the first layer of our approximating scheme, the r layer. Each of the parameters n, e, 
and 5 are present in the notion of solution, Definition 13.31 below, but they are frozen in this section, so we only 
indicate dependence of the approximating sequence on r, the time splitting parameter. We partition the time 
interval [0,T] into — time intervals of length r, where — is assumed to be an even integer. Denoting tj = jr , 
we define the functions h^ et and hf t via 


2 r x 

fcdet(s) = H = 1 - h T st (s). (3.1) 

3=0 

The main result of this section is the following: 

Theorem 3.1. Let (Cl, JF, P, {/3 fc })f =1 ) be a stochastic basis and suppose 1 0 is the filtration gen¬ 

erated by the collection of Brownian motions {P k } k=1 . 
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There exists a sequence {(/3 t , m t )} t >o o/t layer approximations (in the sense of Definitio ns. 3\ below), relative 
to the given stochastic basis, such that for all p > 1 


sup Ef 

T >0 

sup 

r>0 


+ \Pr\7r(LZ) I^det^T| L ^(jyi ) 


< oo 


+ fi 2 P? )l %(L%) 


(3.2) 


< OO. 


Let us proceed to a precise definition of a r layer approximation. To do so, we introduce three elements 
of our approximating scheme: a finite dimensional space where the velocity evolves, a regularization of the 
multiplicative structure of the noise, and an artifical pressure. 

Let {(X„, n n )}()? =1 be a collection of finite dimensional subspaces of [L(fi\ d , together with a sequence of linear 
operators H„ : [L(fi\ d — > X n which satisfy: 

Hypothesis 3.2. The space X n is a spanned by a finite number of compactly supported vector fields in [C^.] d . 
n n : {L(fi\ d —y X n is a linear operator. Let s £ {0,1, 2,3} and 1 < p < oo. For each u £ [W(fi] d 

lim |II n u — u\ w s,d+i = 0. 

n—> oo vv x,o 


The pair (X n ,H n ) can be constructed using a wavelet expansion. For more details on wavelet expansions in 
domains, see [25]. Next, let C+ be the cone of positive functions in C x and ps a standard mollifier. Define the 
operator a k , T ,n,8 ■ C+ x [L].] d ->• X n by 


Vk,T,n,s{p, m) = n n o cr fc (p * rjs(-), [(p A —) — ] * r)s(-), •), 

r p 

where (p, to) are understood to be extended by zero outside of D. 

Finally, the original pressure in the momentum equation will be replaced by an “artifical” one of the form 
P{p) = p 1 + SpP for a sufficiently large power (3. Specifically, we require that 

/3 > max(d, 2y, 4). (3.3) 

Definition 3.3. A pair ( p T ,u T ) is defined to be a t layer approximation to the compressible Stochastic Navier 
Stokes equations HU provided there exists a stochastic basis (Cl T , J>, {T(}J' =0 , P T , {fir}k=i) such that: 

1. The filtration {fF(.}f =0 is generated by the collection {/3 x } k -i- 

2. The pair (p T , u T ) : Cl T x [0, T] —> L(j. x X n is progressively measurable with respect to with P r a.s. 

continuous sample paths. 


3. For all <f> £ C°°(D) and t £ [0,T], the following equality holds P T a.s. 

/ Pr(t)fidx= Po,6<t> + / 2h T det (s)[p T u T ■ \/cj) + ep T X(j)\dxds. 

Jd Jd Jo Jd 

f. For all <f> € X n and t £ [0, T], the following equality holds P r a.s. 

/ p T u T (t) ■ <f>dx = / TOo,5 • cj> + / / 2h T det [p T u T <g> u T — 2pS7u T — A div u T I\ : V0 

JD JD JO JD 

+ f 2h T det [{pl + SpP) divcf) — eVu T Vp T ■ (f>\ dxds 
Jo 


(3.4) 
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n nt n 

+ V'/ / V2h T st p T a k ,T,nApT.PTUr) ■ <t>dxdp k T {s). 

k=1 Jo Jd 


(3.5) 
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5. For all t G [0, T], the following approximate energy identity holds P r a.s. 


id 


p T \u T \ ( t ) + 
z 7—1 


‘Zh'det 




dx 


>0 JD 
n r t 


2p\Xu T \ 2 + X(divu T ) 2 + e(jp) 2 + 5/3$ 2 )|Vp T 

^ J / / Z2h st p T U T • (Jk,T,n,6 )Pt ; PT^"r)dxdf3^. (s) 

fc= lJo 

n rt r 

^ J / / ‘^‘^stP'r I CTk,T,n,5 (/r ; Pt^t) \ dxdt T E n (0) 

, Jo Jfl 


dxds 


fc=i 


sup£' rl ( 0 ) < S 5 ( 0 ) = ^ f 
n J J D 


l m O,<5 | 2 , 1 7 


D L PO,<5 7 - 1 


Po,<5 


da:. 


(3.6) 


In the definition above, we have replaced the initial data (po,mo) by the pair (po,s,mo,s) which satisfy 
Hypothesis 3.4. For eac/i <5 > 0, po .5 €= C°°(D) and 


Vpo ,<5 • n |ai)= 0 0 < d < po,a( 2 ; ) < d 2,3 /or x £ D. (3.7) 

77ie sequence {po,<s}< 5>0 converges strongly to po in the sense that 

Jim |po ,5 - Po\li(d) + \{x £ D \ p o s < poll = 0. (3.8) 

(5-s-O 

The regularized initial momentum (rao,j}i 5 >o ctre defined by the relation 


mo,<5 


mo if Po,s(x) > po(x) 

0 if po,s(x) < po(x). 


(3.9) 


Now we build up to the proof of Theorem 13.11 establishing some preliminary results in Sections urn and 13.21 
below then proving the Theorem in Section 13.31 The essence of the proof is an inductive construction of a r 
layer approximation. On the time interval (0, t], the noise does not contribute to the weak form, so a pathwise 
application of a deterministic result from Section 13.11 will suffice. On the interval (r, 2 t] , the density remains 
frozen and the noise is the sole contribution to the weak form, leading to a simple SPDE for the velocity. 
In Section 13.21 we write this SPDE down and use the regularized multiplicative noise structure to reduce the 
existence to a classical fixed point problem. 


3.1 Machinery from the Deterministic Theory 

Given p £ C+ , define the operator M. [p] : X n —> X* by the relation 

(M[p\u, rj) = / p{x)u{x) ■ p[x)dx 

Jd 

for it, r] £ X n . The proof of the lemma below is left to the reader. 

Lemma 3.5. For each p £ Cff , Aljp] : X n —> X* is an invertible (linear) mapping and 

\ M ~^p\\c{x^x n ) ^ ld _1 |c,- 

Moreover, for each pi, P 2 £ C+ the inverse satisfies the following continuity estimate: 

\M 1 \pi]-M 1 [P2]\ c ^ Xn - ) < \Pi 1 \Cr\P2 ^CxIPl - Pi\l\- 
Let us also introduce the mapping M : C x x X n —> X* by the relation 

{J\f[p, u ], rf) = / [pu ® u — 2pSJ u + (p 7 + Sp 13 — X div u)/] : Vp — eVttV p ■ p dx. 
Jd 
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Proposition 3.6. Let s < T be initial and final times and suppose initial data (pi n ,Um) £ Gff x Cf. are given. 
Then there exists a unique pair (p,u) £ C ((s, T]\C 2 (D) D C + (D)) x C ((s,T]; X„) satisfying the system 


dtp = 2eA p — 2 div(/7u) 


= o 

dn U 


p(s) — Pin 

u{t) = M- l [p{t)\ o (m* n + f*2J\f[u(r),p(r)\dr 
where m* n £ X* is defined for rj £ X n via the relation 

(™*n>v) = / PinUin-V- 
Jd 


in D x (s, T] 
in dD x (s,T] 
in D 
in (s,T] 


(3.10) 


If Ui n £ X n , then u(s) = Ui n . Moreover, the solution map S : C + (D ) x X n —> C ((s, T]; C + (D)) x C {{s,T\,X n ) 
is continuous. 


Proof. The proof uses a straightforward combination of the contraction mapping principle, the L^(Wf ,p ) esti¬ 
mates for the parabolic Neumann problem, the maximum principle, and basic a priori bounds for the system 
(13.101) . A similar result is established in |10] using the Schauder fixed point theorem, though the uniqueness is 
not proven. The details are left to the reader. □ 


3.2 A Classical SPDE Result 

Let (Cl, F, P, {/3fc}fc =1 ) be a stochastic basis such the filtration {J 74 }^ is generated by the collection 

{/3k}k—i- Suppose s < T are two times and p is an T s measurable, C+ valued random variable. An X n valued, 
{.F 4 }^ progressively measurable process u is defined to be a solution to the SPDE 


| d t u = J2k=i &k,T t n,s(p , pu)Pk(t) in (s, T] x D 

|it(s) = u in in D 

provided that for all t, £ [s,T] the following equality(in X n ) holds P a.s. 


« rt 

u(t) =u in + ^2 / Vk,T, n ,s (p, pu(r)) d/3 k (r). 

k —1 ^ s 


(3.11) 


(3.12) 


Proposition 3.7. Given a stochastic basis and a density p as above, there exists a unique solution u £ 
L 2 (Cl; C ([s, T\; X n ) ) to (13.111) in the sense of (13.121) . 

Proof. Define a “random diffusion coefficient” G : X n x!lq X n via 

n 

G{u , (Jj) = ^2 p{u)u). (3.13) 

k= 1 

Note that Hypothesis 13.21 gives L 2 stability of the projections, which may be combined with the continuity 
Hypothesis 11.31 to check that G is lipschitz in u, uniformly in u>. Indeed, for u,v £ X n 

\G(u,w) - G(v,ut)\ Xn < ( [ |[(pM A -)(u-v)] *rj s \ 2 dx)i 

Jd t 

< \p(u>) A —(u — v )\ L 2 < -\ u -v\ L i < \ u - v \ Xn _ 

Here we have used that X n is a finite dimensional space. Hence, the Proposition may be established in the 
classical way via the contraction mapping principle. □ 
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3.3 Proof of Theorem 13.11 


We are now prepared to establish an existence theorem for the lowest level of our scheme. 

Proof. Let (f2,P, {j3k}k=n {F}J=o) be a stochastic basis and assume the filtration is generated by the 

Brownian motions {/3/c}/J = i- We will define the solution inductively. Namely, suppose that (p T ,u T ) have been 
constructed to satisfy the continuity equation (13.41) . the momentum equation (13.51) . and the energy identity (13.61) 
on the time interval [0, t. r 2 j]- To extend the solution to the interval (£ 2 . 7 , £ 2 . 7 + 1 ], apply Proposition 13.61 to find a 
unique pair (p, u) satisfying: 


dtp + 2 div(pu) — 2eA p = 0 

f>£. - n 

dn U 

< p(f2j) = Pt ( t2j ) 

u(t) = M~ l [p{t)} O (p T u T (t 2 j)* + fl. 2Af[u(s),p(s)]ds'j 
u(t 2 j) = Ur{t 2 j) 


in D x (t 2 j,t 2 j+i\ 
in dD x (t 2j ,t 2j+ i] 
in D 

in D x (t 2 j,t 2 j+i\ 
in D. 


(3.14) 


To extend the solution to the interval (£ 2 . 7 + 1 , £ 2 . 7 + 2 ] we appeal to Proposition 13.71 to find a unique pair (p,u) 
satisfying 


dtp — 0 in D x (£ 2 . 7 + 1 ; £ 2 . 7 + 2 ] 

d t u = V2 Y2=i crk,T,n,s(p,pu)/3 k in D x (£ 2 . 7 + 1 , £ 2 . 7 + 2 ]- 


(3.15) 


Using the Ito Formula and the inductive hypothesis, one may check that (ED-flUD continue to hold for 
£ G [t‘ 2 j ,£ 2 . 7 + 2 ]- To prove the uniform bounds, begin with the energy identity (13.61) . One can estimate the 
stochastic integral terms and the Ito correction with the same manipulations as in the formal proof provided in 
Section 12.31 The only additional detail is to note that Hypothesis 13.21 and the uniform boundedness principle 

2 'y 2"y 

imply the projection operators Tl n are bounded (uniformly in n) as linear operators from Lf -1 to Lf -1 . The 
desired measurability, part [2] of Definition 13.31 follows from the continuity of the solution map to the deter¬ 
ministic problem(guaranteed by Proposition 13.61) . together with the fact the that we obtain a stochastically 
strong(measurable with respect to the same filtration as the Brownian motions) solution during each time 
interval where the stochastic forcing evolves. □ 


4 n Layer Existence 


In this section, we apply Theorem 13. II to build the next layer of the approximating scheme, the n layer. Our 
goal is to establish the following: 

Theorem 4.1. There exists a sequence {(/>„, u rl )}^ =:1 of n layer approximations (in the sense of Definition ^- £\ 
below), relative to a collection of stochastic bases {(f 2 „, T n , {.7^}£_ 0 , P n , {$n}k=i)}?^=i> such that for all p > 1 


sup E^ n 

n> 1 


+ \P n \ 1 Lf‘(LZ) 


12 V 

1 £?(»£„) 


< 00 


sup E 1 ^" 

n> 1 



Pp 


|e 5 V(p„ 


■&*Pn)\ 


2 P 

L 2 ALl) 


< OO. 


(4.1) 


Let us introduce the n layer regularization of the multiplicative noise structure. Define an operator ak, n ,s ■ 
L\ x [Lf) d —► X n via the relation 

Vk,n,s{p,m) = n„ o <j k (p * r)s(-),m * r] S (-),-). (4.2) 

Definition 4.2. A pair ( p n ,u n ) is defined to be an n layer approximation to (11.11) provided there exists a 
stochastic basis (f 2 „, JF n , {.7^}|L 0 , P n , {$n}k=i) such that 

1. The pair (p n , p n u n ) ■ x [0,T] —> L% x L%. is an {lFf}JL 0 progressively measurable stochastic process 
with P„ a.s. continuous sample paths. The velocity u n G L 2 x [0, T]; .Hq is an equivalence class of 
{F n }J= 0 progressively measurable Hf valued processes. 
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2. For all (f> £ C°°(D) and all times t £ [0, T] the following equality holds P„ a.s. 

/ p n {t)4>dx = / Po,< 5 ^ + / / [p n M n • V0 + ep„A0]cfa;ds. 
ifl Jb Jo Jd 

3. -For all <j> £ X n and all times t £ [ 0 , T] t/ie following equality holds P n a.s. 

/ PnU n {t) ■ (j)dx = / m 0 l 5 -</) + / / ® - 2pVu n - Adivfi„/] : V0 

JD JD Jo Jd 

+ [(^n + <>/3n) div (j) - eViinVpn ■ (f\ dxds 

Jo Jd 


io Jd 

n ft 


nl /* 

F ^ ^ / / Pn&k,n,8(Pm PnV*n) * f>dxd(3f l (s'). 

k=1 Jo Jd 


(4.3) 


(4.4) 


For all t G [0,T], i/ie following approximate energy identity holds P„ a.s. 


/d 2 7-1 

/•t 


- 1 


P^(i)]dx 


f I [2^|V x w „| 2 + Adivfi 2 + e( 7 Pn 2 + <5/3/3^ 2 )|Vp„| 2 ]dxds 
Jo Jd 
^ /»£ /* 

V / / PnU n - (J k ,n } s(Pn,PnU n )dxd^(s) 

lo Jd 


k=1 

1 

2 


^ ^ / / Pr (^) | &k,n,S (.Pm Priori) | ^^5 -|- ^(O)- 

fc= iJ° Jd 

supF n ( 0 ) < F a ( 0 ) = ^ [ ^ m °’ s 

n ± Jl 


P0,6 7 - 1 


7 

P 0 ,<5 


(4.5) 


For each n fixed we apply Theorem 13.11 to obtain a sequence of r layer approximations {(p T . n , fi Tn )}i->o- 
In Section an we prove a compactness result for this sequence and extract a candidate n layer approximation 
( p n ,u n ) built on a convenient choice of probability space (f2„, F n , P„). In Sectionwe use the compactness 
result to verify (p n ,u n ) is an n layer approximation in the sense of Definition 14.21 


4.1 r —> 0 Compactness Step 

Now we proceed to our compactness step: 

Proposition 4.3. Let W = {/3 fe }fc =1 denote the collection of Brownian motions in the stochastic basis supporting 
our sequence {(p T , 'a T )}r>o- 

There exists a probability space (fl ni F n , P„), a limit point (p ni u n ), and a sequence of “recovery” maps 

{Tt}t> 0 

T t '■ (Lin, J" m Pn) (Ll,iF , P) 

with the properties listed below: 

1. The measure P can be recovered by pushing forward P ra with T t . 


2. The new sequence {(p T ,u T )} T> o defined by ( p T ,u T ) = (p T ,u T ) oT t constitutes a r layer approximation 
relative to the stochastic basis (Cl n , JF n , P n , {iF^}f_ 0 , W T ), where 

W T := W oT t F‘ := a (r t p T , r t u T , r t (p T u T ), r t W T ). 

Moreover, the initial data are recovered in the sense that (p n (0), u n (0)) = (po,s, ■M.~ 1 [po,s]'mo,8)- 
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3. The following uniform bounds hold for all p > 1 


sup E 

r>0 

su P E # 

T> 0 


W~Pt u t\l°°(LI) + \P t \'l°°(l2) + I ^det^r \ (H^ x ) 

+ l^det £2 ^(PT + S 2 Pt )|^(£2) 


< 00 


< oo. 


f. As t —> 0, the following convergences hold P„ a.s. 


Pt ^ Pn 
ILt y Un 

W T ^W n 


in C t (L^)nL^(W^) 
in C t (X n ) 
in [C t {X n )] n . 


(4.6) 


(4.7) 

(4.8) 

(4.9) 


The proof of Proposition 14.31 begins with a tightness lemma. 

Lemma 4.4. The sequence of induced measures {P o (p T ,u T ,W)~ 1 }t>o ore tight on Ct(L f) fl x 

[C t {X n )\ x [C t ] n . 

Proof. Note that it suffices to show the tightness of each component separately. Tightness of P o W~ x is an 
immediate consequence of Arzela-Ascoli and the usual L 2 (Cl; C^ 3 ) bound on each one dimensional Brownian 
motion. Next we will check that 

Jim supP(|n T | Ct(Xn) > M) = 0. (4.10) 

Multiplying and dividing by the density gives the pathwise upper bound 

sup / |m t | 2 < | Pr^L^ sup / p T \u T \ 2 . 
te[0,T]JD ’ te[0,T]JD 

Also, note that if X,Y : Cl —> R are two positive random variables, then 

P(AF > M) < P(A > sfM) + P(Y > VM). (4.11) 


Combining these observations yields 

P( sup f |m t | 2 >M)<P( sup f p T \u T \ 2 > VM) +f{\pf 1 \ L ^> VM). 

te[0,T]JD tE[0,T]JD 


Using the L 2 (Cl) bounds on kinetic energy implied by (13.21) . we can choose M to make the first probability 
small, uniformly in r. 

To treat the second term, recall the splitting scheme from Section [3] defining the evolution of p T . On 
time intervals (t- 2 j ■ fgj+i], Pt solves a parabolic equation with drift u T and remains constant on the intervals 
, £ 2 . 7 + 2 ] • Iteratively apply the maximum principle then use the equivalence of the X n and Cf. norms. This 
controls the second probability from above by 


kjUs° ex P(_^ ft det(*)l div u T (t)\L<x>dt > VM) j 


< 


l^det {t)Ur{t)\x n )dt > C n 1 log(M|p 0 j|JL) . 


Applying a Holder(in time) and the L 2 {X n ) bounds on the velocity implied by (13.21) . we can make this second 
probability uniformly arbitrarily small also. Hence, (14.101) is established. 

We can now bootstrap (I4.I0D and prove the tightness of {P o p T 1 } T> o on C t (Lf) fl I/f (Uj ,/3 ). To this end, 
we use Lemma IB. 71 from the appendix. For simplicity, we will omit dependence of the estimate on the initial 
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density, since it has been smoothed out already. Start by defining the exponent q via the interpolation condition 
i + 2 ( 13 + 1 ) obtain the following estimate: 

\ dt Pr\ L ^ x ~P \P t \l,P(w2' P ) ~ l^det ^ v (PtU t )\l^ x < |&r | C t (X„) \ ^detPr I Lf (Wl' q ) 


< \^r\Ct(X n )\h^ e tP T \ 2 L i3+l \Pr\^i3^ w 2,py 


Applying Cauchy’s inequality, we may close the estimate then interpolate once more to obtain 


\d t p T 


I Pt 


I L^(W 2 ' P ) ~ l^ T lc , t(A„)l^'detPT| i 0+i 1$ \^t\ c t (A' n ) l^detPr I @ ff+i 


< 




~T~ 


< \u 12 


UC t (X n ) 


l^detPr 


'l-t(lZ) 


K e tP P r /2 \ 


< Ifl. 12 


T|C t (A'„)l /l detPr|®co (i g ) |^etPf /2 | 


1-0 

~w 

2d 

lI(lT 1 ) 
1-0 

~W 

L 2 t (wTY 


all 


Note that 6 is defined by the relation ^-j- = 0(1 — |) + (1 — 9). Bootstrapping this estimate once yields for 

r < /3 

|0tVp T |x,r^ + \Y7p T \ L r^ w +r^ < |V div(p T ti T )< | , u T |c t (x„) . (4-12) 

Choosing r large enough to ensure the embedding W*’ r LP is compact, we may conclude from Arzela-Ascoli 
and Aubin-Lions the following set is compact in C t {L&) n Lf(W^ ,/3 ) 

{/ G L*(WW) I |ft/| Lf(ig) + \d t Vf\Lr( L r) + \f\ L!(w; . P) < M}. 

Combining (14.101) together with the uniform estimates (13.21) to control the RHS of (14.121) and (|4.12l) . we obtain 
the desired tightness by Chebyshev for M large enough. Our final step is to show 


lim 

M—> oo 


SUpP([M r ] c l 

T \ t 


1/3 (X„) > M )-0. 


(4.13) 


Note that the brackets indicate we are considering the Holder seminorm, since the uniform norm has al¬ 
ready been handled above. Recalling the operators introduced in Section [3711 define the X* valued processes 
{I?(t)} t > 0 ,{I?(t )} t > 0 via 

I ?(t)= [ 2/idet (rW[u T (r),PT(rj\ dr. 

Jo 

(d?(t),(f>) = / / V2hl t (r)p T (r)a k ,T,n,6(PT(r),pTUT(.r)) ■ cj)dxdj3k(r). 

Jo JD 


>0 JD 

for <j) G X n . For each s < t the momentum equation yields 

u T (t)-u T (s) =M~ 1 [pr(t)\ (4 D W-4 D (s)+l T S (i)-l T S (s)) 

+ {M~ l [p T {t)\ - M-*\pr{s))) o {rn^s + I?(s) + I?(s)) . 

Using Lemma 13.51 and the maximum principle, we obtain the IP a.s. estimate 

K,<slx* “P \d T ]cy 3 ( A'*) "P [d-r 1c)^ 3 (a*) 


(4.14) 


K] c va (Xn) < e°^\ 


+ e c 'n' r l*'rlc t (x„) 


\Pr\Cl /3 (Ll) {\d?\c t (X*) + Urlctix-)) 


(4.15) 


In view of the estimates for the density above, this reduces the problem to controlling the probability that Zf 
and 1^ have a large Holder norm. To estimate \I®\ c i/ 3 ^ x ,y note first that for p G L], and u G X n 

|A/'(p,u)|a-* < Cn (\p\Ll\u\x n + \u\x n + \p\ly +S\pf L p + \u\xjp\ w + i) • 
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Applying Holder’s inequality in time yields for all s < t 


J |^(p T (r),u T (r))|x* T> [\Pr\L^(Ll)\u T \c t (X n ) + I^H~(L2) + l^oo^j 

+ (f-s) 1 ~? \u T \ Ct (x n )\p T \ L P( W ^) ■ 

Certainly 1 — ^ Hence, we may combine (Id. 101) . (14.121) . and the uniform bounds (13.21) to obtain 

lim supP > m') = 0. 


M—too T 


(4.16) 


To estimate |/,f Li/3 , v ,., fix a </> € X n . Apply the BDG inequality, the boundedness of the projections, and the 
summability Hypotheses 11.41 for the noise coefficients to obtain for all p > 2 

r t 


E p 


< E p 


a JD 
t 


V2hl t (r)p T (r)a k , T ,n,s{pT{r),p T u T (r)) ■ (j)dxdj3 k {r)\ p 


< Wk\ p 


( / p T (r)a k ,T,n,s(pT(r), p T u T (r )) • 4>dx) 2 dr\ * 
L Js JD 

E p ' 


C " 1 


'[I [ \Pr(r)\h\ p/2 } <(t- S ) p/2 supE p [|p r | 
J s r 


V 1 


This yields for all s < i, p > 2 and <p G X n 


supE p [|(/f,^)|^ s , P ] < supE F [|/5 T 


r>0 


r> 0 


T r (L 2 )J- 


Choose s,p such that Wf’ p C\^. Since X n is a finite dimensional, the uniform bounds (13.21) imply 

i} lS M'\x*j > M )= o- ( 4 - 17 ) 

Starting with the identity (14.151) and using (14.161) . (14.171) and some elementary estimates similar to (14.111) give 
the tightness of the laws {Po u“ 1 } T>0 on C t (X n ) by Arzela-Ascoli. □ 

Next we apply the tightness result above together with a version of the Skorohod Theorem lA.12l to complete 
our compactness step. 

Proof of Proposition El : Let us define the sequence of random variables {X r } T> o via the relation 

X T = {p T ,u T ,W). 


These random variables induce a tight sequence of laws on the metric space 

E = C t {L^ x )nL p t {W 1 J) x [C t (X n )\ x [C t (X n )] n 

by Lemma [4.41 In view of Remark IA.131 (E,t) is a Jakubowski space and we may apply Theorem IA.12I to 
obtain a new probability space (Q n ,J r n . P n ), a sequence of recovery maps {T t } t> q and a limiting random 
variable X n = (p n ,u n , W n ) such that parts [l] and |4] of the Proposition hold. To check the uniform bounds (|4.6|) 
on ( p T , u T ), simply use the recovery maps together with the bounds on the original probability space (13.21) . This 
is rigorous because the following functional is is continuous from (E,t) to (R, S(R)) and hence measurable. 


(p,u,W) -»• \VP u \ 2 L P r ( L i) +\pVl T {lZ) 
+ l^det e 2 V(p 2 + P*)\l2(lz)- 


\&P\ 


ftp 


+ l^det M l 


2 P 

L 2 (Wp 2 ) 


(4.18) 


To check part 2, note that P(T T (f2 n )) = 1, so the recovery maps allow us to preserve the continuity equation. 
Recovering the momentum equation and checking adaptedness with respect to the new filtration requires a 
small argument regarding the recovery of stochastic integrals. This is not used directly in the proof below so 
we omit it, for more details see n. □ 
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4.2 r — ► 0 Identification Step 

Recall that r t denotes a restriction operator. Consider r* on the spaces dictated by Proposition 031 and define 
a filtration {7™}|L 0 by 

Pf = <j(r t Pn,rtUn,rt(pnU n ),r t W n ). (4.19) 

Lemma 4.5. The pair ( p n ,u n ) satisfies the parabolic equation, (14.31) of Definition \f.^ 


det ~^ 2 1 

□ 


Proof. By part [2] of Proposition 14.31 the parabolic equation (13.41) is satisfied by (p T ,u T ) on Cl n a.s. with 
respect to 1P„. Applying the pointwise convergences (l4.7H - d4.8l) along with the weak convergence of h T - v 1 
we easily pass the limit in the weak form. 

Lemma 4.6. The pair (p n ,u n ) satisfies the projected momentum equation, (14.41) of Definition m 
Proof. Given cj> G X n , we define the adapted, continuous stochastic process {M t n }^ 0 via 


D 


b 


/ p n u n (t) ■ <j)dx — / mg ■ (j>dx — / / [p n u n ® u n — 2pS7u n — A div u n ] : V<^> dxds 


o Jd 


[ [ [{Pn + SPn) di v ~ fVu n Vp„ ■ f>] dxds. 

Jo JD 


In an analogous manner, we define the processes {M[{(/>)}JL 0 and {Mf (<f>)}-jL 0 via 

Mf((f>)= / p T u T {t) ■ <j)dx - / TOq • (fdx — / / 2h Aet [p T u T ® u T — 2pS7u T — A div u T ] : Vcfdxds 


o Jd 


+ f ( 2h T det [{pf + SpP) div (j) + 2eVfi r V/3 T • </>] dxds. 

Jo Jd 

Mf (fi i) = / p T u T [t) ■ (f>dx — / TOq ■ fidx — / / 2h det [p T u T ® u T — 2pS7u T — Adivu r ] : V0 

Jd Jd Jo Jd 

+ f f h Aet [(pf + 5pP) div + 2 eV u T Vp T • <f\ dxds. 

Jo JD 


io Jd 

Our plan is to check the criterion laid forth in Lemma 1A. 151 in order to identify 

n ft 

1 0 JD 


_ n ft f 

Mjf((f>) = V' / / PnVk,n,s(p n ,P n U n ) ■ (t>dxdj3*(s). 

k=! J o J D 


This implies the momentum equation (14.41) holds. Let us fix in advance two arbitrary times s < t and a 
continuous functional 7 

7:C([0 ,s];L£) x C([0, «];*„) x C([0, s]; if) x C([0, s]; X n ) -> R 

which will be used repeatedly below. We will repeately use the fact that in order to verify a process {N t }]L 0 is 
a martingale on (fl n , T n ,F n ), it suffices to verify 


E r 


l{r s p n ,r s u n ,r s (p n u n ),r s W n )(N t - N s ) 


= 0. 


We start by using the Levy Characterization to verify {j3*(t)}J =0 is an (71* }^1 0 Brownian Motion. Applying 
the pointwise convergences (I4.7D - (14.9I) together with the uniform bounds (14.61) we find that 


E f 


'y(r s p n ,r s u n ,r s (p n u n ),r s W n )(f3*(t) - $%(s)) 


= lim E r " 

r —>-0 


7 (r s p T ,r s u T ,r s (p T u T ),r s W T )(/3f(t) - /3^(s)) 
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Using Part [T] of Proposition 14.31 with a change of variables, then the martingale property of 0 k , we deduce 

7 (r s p T , r s u T , r s (p T u T ), r a W T ) (/3*(t) - /3*(s)) 
l{r s p T , r s ii T , r s (p T u T ), r s W) (0 k (t) - 0 k {s)) I = 0 . 


lim E p " 

r->0 


= lim E F 

T-> 0 


Similarly, one verifies 


E p 


l(r s p n , r s u n , r s p n u n , r s W n ) (/3^(t) 2 - /3^(s) 2 - i + s) 


= 0 . 


Next we check that {Mfi ((/))}f =0 is an {.F "}£_ 0 martingale with quadratic variation 

\ 2 


^ ^ Pn&k t T,n,5(,Pn: Pndtn) ‘ 4*^ ds. 


Recall that /ij et —*■ i weakly. Hence, by using (14.7D - fl4.9H together with (14.61) ; followed by (13.51) of Definition [I] 


we obtain: 


E p 


7 (r s Pn,r s u n ,r s (p n u n ),r s W n )(M?(<j)) - M" (</>)) 


= lim E p 

T->0 


7 [r s p T ,r s u T ,r s (p T u T ),r s W T ){Mf((f)) - 


lim E p 

r—>-0 


7(' r s p T ,r s u T ,r s {p T u T ),r s W T )(M[(4 1) - 


= 0 


In the remaining analysis we will suppress the arguments of 7 . Similarly, using the same facts as above along 
with the weak convergence hf t —> ^ we find 


= lim E p 

T— >-0 


= lim E p 

r—>0 


E Fn (M"(^)) - (M g n (^))“ - X] y (J PnVk,T,n,s(p n , PnUn) ■ (frdx) 2 dr 

1 7 2 hl t (r)(J^p T a k<Ttnt s(p T ,p T u T )-(j)dx) 2 dr 

7 i^J^t (0)) 2 - { M l{<t>)) 2 - HZ [ 2h st( r )( [^PrCTk,r,n t s(Pr,PrUr) ■ (fdx) 2 dr 
In essentially the same way, we check 


fc=1 1 


id 


= 0 . 


E p 


1 ((/)) 0n(t) — M"(</>)/3*(s) — J p n a k ,r,n, 5 (p n ,u n ) ■ (j)dxdr 

Using Lemma [A. 151 we conclude. 

Lemma 4.7. The pair ( p n ,u n ) satisfies the energy identity, (14.51) from Definition ^ .2\ 
Proof. Define the {7’*}^1 0 continuous, adapted stochastic process {Et}t=o y i a 

K= I []-p n \u n \ 2 (t) + ^—pZ(t) + -^-p^t)]dx 

Jo 2 7-1 P -1 


= 0 . 


□ 


+ [ [ [2p\Vu n \ 2 + X(dwu n ) 2 + e{pl 2 + <5/5£ 2 )|V i o n | 2 ] dxds 
Jo Jd 

n rt r 

^ ^ / / Pt (^) |^/c,r,n,i5 (Pr) Pr Hr ) I dxds U n (0) 

do JD 
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We may now implement the same method as in Lemma 14.61 in order to identify o as 

n i*t r 

^ ^ / / Pn'U’n ' ^fc,T,n,(S (Pr j Pt') dxd(3 r (s). 

fc=lV° 

One can now perform the same manipulations as in the formal estimates in order to deduce the uniform 
bounds. □ 

Proof of Theorem \4-l\ For each n > 1, we obtain an n layer approximation ( p n ,Un ) using our compactness 
step, Proposition 14.31 Indeed, we can check each part of Definition 14.21 as follows: Part [1] follows immediately 
from the definition of the filtration {J-f}J = 0 , the continuity and momentum equations are satisfied in view of 
Lemmas 14.51 and 14.61 while the approximate energy identity 14.71 follows from Lemma 14.71 To obtain the uniform 
bounds 14.11 use the energy identity along with the method described in Section 12.31 □ 


5 e Layer Existence 


This section is devoted to the e layer existence theory; sending n —> oo our goal is to prove: 

Theorem 5.1. There exists a sequence {(p e ,u e )} e> o of e layer approximations (in the sense of Definition tS7^\ 
below), relative to a collection of stochastic bases {(D e , F e , {E(}(L 0 , P e , {/3 £ }kLi)}e>o, such that for allp> 1 


supE^ 

£>0 

sup 

£>0 


I y/Ts 


12 V 


e 'L?(LZ) 

L i Pp 


\u, 


12 V 


Wrw 

L»(i|) + V ^ 2 + S *P‘ )l %{Ll) 


1 -Sr^2 P 


< OO 


< oo. 


(5.1) 


For each fc, introduce the operator a k ,s : L) x [Lf] d C ° via the prescription 


<Jk,s(p,m) = <7 k {p * r75(-), m * Vs(-), ■)• 

Definition 5.2. A pair ( p e ,u e ) is an e layer approximation to (11.11) provided there exists a stochastic basis 
0 . p ^{^} fc > 1 ) such that 

1. The pair (p e , p e u e ) : x [0,T] —► [Lf x Lj ¥I ] w is an 0 progressively measurable stochastic proces 

with P e a.s. continuous sample paths. The velocity u e £ L 2 x [0, T]; belongs to the equivalence 

classes of {E(}JL 0 progressively measurable Hq x valued processes. 


2. For all f> £ C°°(D) and all t £ [0, T], the following equality holds P e a.s. 

/ p £ (t)<j>dx = / po,6 + / / [PeUe ■ V0 + ep e Acf]dxds. (5.2) 

J d Jd Jo Jd 

3. For all (j) £ [Cf°(D)] d and all t £ [0,T], the following equality holds P e a.s. 

/ p e ii e (t) ■ <f>dx = / mo,s-4>+ / / [p e u e 0 u e — 2pS7u e — A div u e I] : S7tfdxds 

Jd Jd ’ Jo Jd 

+ [ f [(pf + SpP) div - eVufVpe ■ <j>\ dxds ( 5 . 3 ) 

Jo Jd 

00 pt p 

+ V / / p e CTk,s(Pe,PeUe) ' 4>dxd^ (s). 

Jo Jd 


For each e > 0 fixed, we can apply Theorem 14.II to obtain a sequence {(/S ra>e , u n:e )}(f =1 of n layer approxi¬ 
mations satisfying the uniform bounds (14.11) . In Section 15.11 we switch probability spaces and use the recovery 
maps to define a new sequence {{p n ,e,u nte )}^ =1 and obtain compactness. We extract a limit point ( p e ,u e ) then 
verify our candidate is an e layer approximation in Section 15.21 


20 






5.1 n —> oo Compactness Step 


Next we establish the following compactness result: 

Proposition 5.3. There exists a probability space (f2 c , P e ), along with a sequence of recovery maps {T„}^L 1 
and limit points (^p e , u e , pplu e , pi + 

T n ■ (Cl e ,P £,P e ) <Sl n ,P n,P n ) 

such that the following hold: 

1. For each n, the measure P n may be recovered from P e by pushing forward T n . 


2. The new sequence {(p n , Un)}^^ defined by (p n ,u n ) = (p n ,Un) °T n constitutes an n layer approximation 
relative to the stochastic basis (fl e ,P e , P e , {Ffyf—Q, W n ), where 

hPn •— O T n J~ n .— O {r t p n , T't'Oni Tt ( Pn'On) ‘Tt lErc) 

3. The following uniform bounds hold for all p > 1 


supE** 

n 

sup E^ 6 

n 


\VPn u n\Lf>( L %) + \Pn\7 rW ) + l Un l4( ff o„ 


2 P 


< OO 


1^ + l £2 ^{Pn + Pn J\l*(L%) 


12 p 


< OO 


4-. The following convergences hold P e a.s. 

Pn Pe 
u n -»• u e 

Pn’O'n t f>eUe 

w n ->w e 

5. The following additional convergences hold 

\JPnU n \Z~feUe 
U n U e 
Pn Pe 


in C t ([LP\ w ) n L 2 t (L 2 x ) 
in [L 2 (Hl x )} w 

2/3 

in C t ([Lf +1 ] w ) 
in [C t }°° 


m L p w ,{n t ;L?(Ll)) 
i n L 2 {Hq x )) 

in L p w , (SX; L?{L%)) n L p w ((l t - L 2 (W l x 2 )) 


(5.4) 


(5.5) 

(5.6) 

(5.7) 

(5.8) 


(5.9) 

(5.10) 

(5.11) 


2/3 

The first step is another tightness lemma. Enumerate a dense subset {c/)k}kLi of Lf +1 . Define the sequence 
of random variables {(X n ,Y n )} n > i via the prescription 


X n — {^Pni U ni p n Uni {P n }k=l'j 
Y n = {(pnUn, 4>k)L^. }fc=l- 


Our convention is that given a topological vector space G , a finite sequence {xfc})J =1 is viewed as an element of 
G°° where Xj = 0 for j > n. These random variables induce measures on the following topological spaces 


E = C t ([L%\ w ) O L\ x x [L 2 (Hl x )\ w x 
F = [C t }°° . 


T OO 


2/3 

(LP 1 ) 


x [G t ] c 
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Lemma 5.4. The sequence of induced measures {P„ o }X n ,Y n ) 1 } n >\ is tight on E x F. 

Proof. It suffices to consider each component of X n separately. The tightness of {P n o (u ni /5 n u n ) -1 }„>i 
follows immediately from the bounds (14.11) and Banach Alaoglu. To treat the collection of SBM, note 


sup E Pn [|/3^| 2 i] < oo. 

k,n: k<n Cfr 


For each M > 0, the set below is compact in [Ct]°° by Arzela-Ascoli and Tychnoff. 


(5.12) 


l[{f € C t | |/| i < M2 3 } (5.13) 

3 =1 * 


Choosing M > 0 appropriately and summing a geometric series gives the desired tightness of {P„ o Wy 1 } n >i. 
Recall that /5„(0) = pa,8 by Part [5] of Proposition 14.31 Since /3 > d we may choose a 9 £ (0,1) and define a 
q > 2 by the relation | = 9{^ + \ + (1 — 9){j^ + |). Maximal regularity results for parabolic equations 

and interpolation give the P„ a.s. inequality 


\dtPn\L q t (w~ 3 ' q ) + e \P™\Ll(wt q ) ~ \poAl%, + 

+ 1 + |/nlln| 2*/3 |Pn^n| 20 

LFiLlX 1 ) 


< 


1 + |/5, 


15 ( 1 + 0 ) 1 


\i-e 


\ I \/Pni 


,i-e 


(5.14) 


Hence, the LHS is uniformly controlled in L 2 (fl n ) in view of the uniform bounds (14.11) and Holder(in w). 
Corollary I A. 1 71 and Aubin-Lions imply that for each M > 0, the following set is compact in C t {[L^.] w ) D L 2 X . 


{feL?(LP)nL 2 \\d t f\ 


I Lt{W- 1 ' q ) 


I/I 


L q t {wT q ) 


+ 1/1 




< Vm} 


Using the uniform bounds on {p„}^ =1 in L 2 (H„; L^°(Lf)) gives the tightness of {P n o / 5“ 1 }^T 1 for an appropriate 
choice of M >0. To address the sequence P„ o Yff 1 , let s < t and <f> £ X n be arbitrary and use the momentum 
equation (14.41) to decompose ( p n u n (t ) — p n {s)u n (s), 4 >)lI into three terms: the stochastic integrals, the energy 
correction, and the rest. To estimate the stochastic integrals, we use the BDG inequality together with the 
stability of the projection operators, Hypothesis 13.21 via 


E 


L fc=l 


< E J 


L k =l 

OO 


( 5 ^/ / Pn<?k,n,s{pn,PnUn) - (/>dxd(3n) P 

: Js J D 

oo ft r 

(^2 j (/ PnCrk,n,s(Pm PnU n ) ■ 4>dx) 2 dr)‘ 

LUL~ m )' “ (/ \Mr)\lz dr ) 


(5.15) 


?==i 

if UIP 


To estimate the energy correction, we use (15.141) to obtain the following inequality P„ a.s. 


n Vu n Xpn ■ (fdxds | < e\(j)\ L c™ / |Vu„(r )| L 2 \S7p n {r)\ L 2 dr 

3 4s 


< e(t-s) 2 jlVpnli^. 
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To treat the remaining terms, Holder and Sobolev yield the P„ a.s. inequality 


/ / [PnU n ®«n- 

Js JD 

2/iVu + (/% + 5f>n — divu„)/] : V (j>dxdr 


(Wx’ 2 ) (1 

+ C(t-s)* IWIloo 

J P n \l?(L2) + ■ 


Combining (I5.15I) - (I5.16I) and using (Id.II) together with (15.141) we obtain for all k and p > 1 
SUp [\(p n Un(t) ~ PnUn{s),<j>k)L* H < C\(j>k |^i \t - s| p( W } . 

n>k 


(5.16) 


(5.17) 


Combining this observation with the Sobolev embedding theorem for fractional sobolev spaces(in time), for any 
a < A — there exists a p such that 

z q 


sup E p " [\{p n U n Ak) C^\ V ] % I&I& 1 . 

n>k x 

For each M > 0, define the set Km via 


(5.18) 


j =1 

In view of Arzela-Ascoli and Tychonoff, Km is a compact set, and Chebyshev yields 


(Y n ? 


K 


M 


< 




fc =1 


PnU n Ak)cf > M p 2 P \(j> k \ C l 


n 



< AT 1 . 


This implies the tightness of the sequence {P n o Y n 1 }„>i. □ 

Now we can complete the proof of our compactness step. 

Proof of Provosition \5.?h : In view of Remark I A. 1 .11 E x F is a Jakubowski space. Thus, we may apply 

Theorem IA.12l in order to obtain a sequence of maps {T„}^L 1 

T n : (Cl e ,fF E,Pe) —> {Cl n ,fF n,Pri) 

and a limiting random variable (A e , Y e ) = (p e , m e , u e , W e , | ). Moreover, the properties listed in Theorem 

IA.12I imply directly Part [1] of the Proposition and guarantee that 


{pn ; Pn^m bFn: fyk)} k—l) — {Pni Priori : 'Wm { {Priori: ^fc)}fc=l) ° 

= (pn°T„, «„oT„, Pn°T n U n of n , {^of n }^ =1 ,{(p n tt n of n ,(/lfc)}^ 1 j (5.19) 

= (X n ,Y n )°f n ^(X e ,Y e )- 


The limit is understood to hold P e almost surely in each of the topologies where the tightness was proven. 
In particular, we obtain the pointwise convergences (15.51) . (15.61) . and (15.81) . To obtain & we need a few 

2/3 

more arguments. By (15.191) . p n u n converges to fh e in [L^°{LS +1 )] W * P c almost surely. Moreover, for each k,n 
satisfying k < n, ( p n u n , 4>k) converges to Y^ in Ct- Hence, for each fixed w in a set of full P e measure, we may 

2 f< 

appeal to appendix lemma IA.16l in order to obtain the convergence of p n u n {ui ) to m e (w) in Ct([L£ +1 ] w ). Indeed, 
the first remark yields the necessary uniform bounds and the second gives the equicontinuity. It is important 
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to remark that since we identified the limit point rh t in advance, there is no need to pass to a subsequence 
and we do not run into any trouble with picking different subsequences for different w. Note that (15.51) implies 
that p n p e strongly in Lf(W~ 1,2 ), which may be combined with (15.61) to identify rh e = p e u t in the sense of 
distributions, P e almost surely, yielding (15.71) . 

It may be checked with a regularization argument that the energy functional defined in (14.181) remains measurable 
with respect to the new topology introduced in this section. Hence, we may recover the uniform bounds (15.41) 
from the prior probability space. Finally, we may use these bounds along with Banacli-Alaogolu theorem to 
obtain &(EIID ‘ □ 


5.2 n —> oo Identification Step 

Consider the restriction operator on the spaces dictated by Proposition 15.31 and define a filtration by 

Tl = o{r t p e ,rtUe,rt(peUe),rtW e ) (5.20) 

Lemma 5.5. The pair (p e ,u e ) satisfies the parabolic equation, (15.21) of Definition [572[ Moreover, we have the 
following convergence upgrade: for all p> 1 

lim ¥F e [\p n -p, 


e,P L?(wfi 2 ) 


= 0 


(5.21) 


Proof. In view of Part [2] of Proposition 15.31 (p n ,u n ) satisfy the same parabolic equation on the new probability 
space f l e almost surely with respect to P e . Appealing to the pointwise convergences (ES-dSU), we may pass to 
the limit in the weak form P e a.s. and obtain the same equation for (p e ,u e ). To prove the convergence upgrade, 
begin by appealing to Lemma iRbl and obtain the following energy identity for all t € [0,T], P e a.s. 


/ p 2 n (t)dx + e / / \Wp n \ 2 dxdt = / pl s dx— / / div u n p 2 n dxdt. 

Jd J o Jd Jd ’ Jo JD 

[ p 2 (t)dx + e f j \Vp e \ 2 dxdt= f Po 6 ~ I f divu e p 2 dxdt. 

JD Jo JD JD ' Jo JD 


(5.22) 

(5.23) 


/0 JD 

Using again (I5.5f) - (|5.6[) . we can pass limits on the RHS of (15.221) and conclude from the LHS of (15.231) that 
a.s. 


lim e 

n—^oo 


0 JD 


|V p n \ 2 dxdt = e / / \S7p e \ 2 dxdt. 

Jo J D 


(5.24) 


Combining this observation with the pointwise convergence (|5.5I) and the uniform bounds (14.11) , one obtains 

}fi^\P n \ L p(p. e -L 2 t (Wfi 2 )) ~ ^^Lp[n e -L 2 t {wfi 2 ))' 

Hence, we may upgrade the weak convergence (15.71) and obtain (15.211) as desired. □ 

Lemma 5.6. The pair (p e ,u e ) satisfies the energy corrected momentum equation 1 5. Si from Definition } 5. 2\ 


Proof. We follow the same general strategy as in Lemma 14.61 Indeed, first note that the exact same proof 
works in order to check that a collection of 0 independent brownian motions. Next, for each 

f> G U£° = 1 X n we introduce a continuous {T^}J =0 adapted stochastic process {M t e (0)}^ o defined by the relation 

r t 


Ml((j)) = / PeUeit) ‘ <\>dx ~ / 7710 ,( 5 - 0 - 


ID 


ID 


[Jv. 


— 2 /rV'Ue — Adivu e /] : Vfidxds 


o Jd 


\(pl + ) div 4> — eVu e Vp e ■ </>] dxds. 


In view of the pointwise convergences (15.51) - (15. 61) . the following limits hold P € a.s. for all t € [0, T] 

lim / p n u n (t) ■ 4>dx = / p e u e (t) ■ fidx (5.25) 

n ~^°° Jd Jd 

lim / [ \2pS7u n + A div u n I] : 'S/ifdxds = [ [ [2/rVu e + A div u e ] / : 'Vfidxds. (5.26) 

"-’Wo Jd Jo Jd 
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Noting the compact embedding L ^ ^ W x 1,2 , we may upgrade (15.71) with Lemma lB.ll and obtain P c a.s. 
p n u n p e u e in L 2 (W~ 1,2 ). Combining with (15.61) we have a weak/strong pairing and obtain P e a.s. for all 

te [o,r] 


lim 

n—> oo 


p e u e <S> u e : 'S/cfdxdt. 


(5.27) 


p n u n ® u n : V(j)dxdt = 

/o Jd Jo Jd 

Combining (15.21|) . the strong convergence upgrade for the density, with the weak convergence of the velocity 
(15.61) . yields(along a subsequence) P e a.s. for all t G [0,T] 


lim 

n—>• oo 


e'Vu n VPn • <t>dxds = 


0 JD 


eVu e S7p e ■ 4>dxds. 


(5.28) 


0 JD 


Recalling the interpolation argument in the proof of the r —> 0 tightness Lemma 14.41 we may use the uniform 
bounds S3 to obtain further: for all p > 2 


supE F '|p n |^ +1(L g +1) < oo. 

Combining this observation with (15.51) gives P e a.s. for all t £ [0,T] 


lim [ I (pi + 6pn) div cfdxds = [ [ (pi + 5p (!) div cfdxds. 
i-s-oojo Jd Jo Jd 


(5.29) 


(5.30) 


These remarks allow us to proceed as in Lemma T4.61 and conclude that {M t £ (</>)}^b 0 is an {Tl}f =0 martingale. 
Indeed, a dominated convergence argument reduces the problem to checking that for all k > 1, P £ a.s. for all 
t G [0,T] we have the convergence 


lim 

n—>-oo 



Pn(t)cr k ,n,s(Pn(t), PnU n (t)) ■ <t>dx 



Pn(t)ak,s(Pn(t), PnU n (t)) ■ (j)dx. 


(5.31) 


2/3 (3 

To check this, note first that a k ,5 : L& x [Lj^ ] d —> L is a compact operator. Since o k n ^ = H n o a k ,s we 

may use the stability of the projections in L]^ 1 and (15.51) . (15.61) in order to conclude that P e a.s. for all t € [0, T\ 

p 

we have ak,n,s(pn(s), p n u n (s)) —> crk,s(p(s),pu(s)) in Lf -1 . Hence we have a weak/strong pairing and we may 
conclude. 

With these additional convergences at hand, it is straightforward to implement the method in Lemma 14.61 
and identify 


OO n 

Mt(<t ) ) = / a k,s(Pe, Petit) ■ <l>dxdj3e(s)- 

k= 1 J ° J ° 


Proceeding by a density argument, using Hypothesis (13.21) and appropriately redefining the stochastic integrals 
on sets of measure zero if necessary, we obtain a single P e full measure set where the energy corrected momentum 


equation holds for all <j) € [ C°°(D )] and t € [0,T]. 


□ 


5.3 Conclusion of the Proof 


Proof of Theorem \5.U For each e > 0, we obtain an e layer approximation ( p e ,u e ) using our compactness step, 
Proposition 15.31 together with Lemmas 15.51 and 15.61 In view of (15.91) - (15.1 ID and the lower semicontinuity of the 
relevant norms, we obtain for each e > 0 


E p 


lv^ 


12 p 

e 'L?(LZ) 




I 


< lim inf E r ' 

n—>oo 


\V~Pn u n\L°o( L 2) 


Kllrw) 


L?(l£) ^ 1 e| Lf(Wj[’ 2 ) 
1 6*p< " 


LF(L%) 


n 'LZ(Wf 2 ) 


Moreover, in view of (15.5I) - (I5.6I) we may identify y/~plu e = y/plu e in the sense of distributions P e almost surely. 
Appealing to (14.11) . we may maximize the inequality in e and obtain the uniform bounds (15.11) . □ 
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6 S Layer Existence 


In this section, we build our final approximating scheme; the <5 layer. Sending e —> 0, our plan is to prove: 

Theorem 6.1. There exists a sequence {(ps, W( 5 )}<s>o of 5 layer approximations (in the sense of Definition \6.2\ 
below), relative to a collection of stochastic bases {(( 2 , 5 , JFg, {Tg}^L 0 , P 5 , {$$ }fcLi)}< 5 > 0 ; such that for allp> 1 


supE n 

< 5>0 


W%' 


su s \% 


(LI) 


I M] 


\S p Ps 




+ \us\ 2p 




< 00 . 


( 6 . 1 ) 


Let us proceed to a precise definition a of 6 layer approximation. 

Definition 6.2. A pair (ps, us) is defined to be a 5 layer approximation to (11.11) provided there exists a stochastic 
basis (Qs,^5,{^s}J=o^s,{Ps}T=i) such that 

1. The pair (ps, psus) : ( 2,5 x [0,T] —>■ [L% x is an -q progressively measurable stochastic process 

with a.s. continuous sample paths. The velocity us £ L 2 ^ 12,5 x [0,T]; Hq belongs to the equivalence 
classes of {Jg}f = 0 progressively measurable x valued processes. 


2. For all f> £ C£°(D) and for all t £ [0, T] the following equality holds P 5 a.s. 

/ ps(t)(j)dx = / po,s-f>+ / / psus ■ Vfdxds. 

JD JD J 0 JD 

3. For all 4> £ [Cf°(D)] d and for all t £ [O,! 1 ], the following equality holds Pj a.s. 

/ psus(t) ■ 4>dx = / mo,j • <p + / / [psus ® us — 2pVus — A div us I] : Vtpdxds 
J D JD J 0 J D 

pt p 00 pt p 

+ (p's + tips) div 4>dxds + Y, / / ps<Jk,s(ps,Psus) ■ (/>dxdj3g(s). 

J 0 J D 


( 6 . 2 ) 


(6.3) 


For each 6 fixed, we can apply Theorem 15.II to obtain a sequence {(p e ^, M e /)} e >o of e layer approximations 
satisfying the uniform bounds (15.11) . In the next section, we prove Proposition 16.31 which allows us to switch 
probability spaces and use the recovery maps to define a new sequence {(p £ ,s, Ue,, 5 )}e>o and obtain compactness. 
However, in order prove tightness of the pressure sequence, we need further estimates on the moments of the 
density. This is proved in Proposition ^. 41 Next we extract a limit point (p e ,u e ) and work to verify our candidate 
is an e layer approximation. This involves both a preliminary identification step, Lemmas 16.61 and 16.71 and an 
elaborate procedure(modelled on the work of Lions [2T|) for proving the strong convergence of the density in 
Lemmas 16.81 and 16.91 


6.1 e —> 0 Compactness Step 

The main goal of this subsection is to prove the following compactness result: 

Proposition 6.3. There exists a new probability space (fls, IFs, Pa), along with a sequence of recovery maps 
{T e } e> o and limit points (^ps, us, y/psus, Ps + $Ps > PS log PS, ps div us^j 

f e :(£ls,fs,h)^(n e ,f £ ,F e ) 

such that the following hold: 

1. The original probability measure P e may be recovered from by pushing forward T e . 
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2. The new sequence {(p e , w e )}e>o defined by ( p e ,u e ) = (p e ,uf) o T e constitutes an e layer approximation 
relative to the stochastic basis (£lg, {J-f}f =0 , W e ), where 

W e := W e o T e T\ \= a (r t p e ,r t u e ,rt(peUe),rtW e ). 


3. The following uniform bounds hold for all p > 1 


supE P5 

e>0 




|2 p 

' e 'L?(Ll) 


17 V 


+ l^lLf (LJ) 


I ^ Pe 


L?(Li) 


12 p 

e 'Ll(Wf’ 2 ) 


< oo. 


f. The following convergences hold a.s. 

Pe ->■ ps 
u e us 

p e u e psus 

pf + 6pP p] + <5/5f 

W £ Wi 

5. TTie following additional convergences hold 

sfplUt -> Vfsus 
Pe PS 
U e —> US 

p e div u t —> ps div us 
Pe log Pe PS log PS 


in C t ([L^\ w ) 
in [L*(W^)\ w 

2/3 

in C t ([Lg +1 ] w ) 


*" Ll(^LT{Ll)) 

m ir wt ((ls-,L?(lg)) 

in LV w (Cls-,L 2 t (W^)) 

2/3 

m Ll{^L 2 t {Lt 2 )) 
m L p w ,(Qs;L?(Ll)). 


(6.4) 


(6.5) 

( 6 . 6 ) 

(6.7) 

( 6 . 8 ) 

(6.9) 


( 6 . 10 ) 

( 6 . 11 ) 

( 6 . 12 ) 

(6.13) 

(6.14) 

(6.15) 


To prove the tightness, we start by proving the following integrability gains: 


Proposition 6.4. The following estimate holds uniformly in e > 0 for all p > 1 


sup E Pe 
£>0 



+ 5pP +1 dxdt\ p 


< oo. 


(6.16) 


Proof. For regular domains D , one can define a sort of “inverse divergence”, known as the Bogovski operator 
B. The properties of B are recalled in appendix lemma IB31 Define the following “random test function” 


<p e = B[p e - j p £ dx\. 


(6.17) 


The parabolic equation and the dirichlet boundary condition for the velocity yields the IP e a.s. equality 

dtfie = eVpe - B[div(p e u e )]. (6.18) 


Since the weak form of the momentum equation is stated in terms of deterministic test functions, “testing” <p e 
requires an appeal to a version of the Ito product rule. The equality below can be justified with a somewhat 
lengthy, but straightforward regularization argument (which we omit) in the spirit of [18j or [5]. For all times 
t € [0, T) we have lP e a.s. 


/ peUeit) ■ ip e (t)dx = / mo ,5 • (0)dx + / / p e u e ■ d t f>e + [peU e ® u e — 2p\7u e ] : V<p e dxds 

Jd Jd Jo Jd 

+ / / (pf + SpP — A div u t )I : V<p e — eVu e Vp e • (p t dxds (6.19) 

Jo Jd 


io Jd 

oo r t 


k=l 


+ / Pe&k,s(Pe,PeUe) ■ Pedxd(3 1 f(s). 


0 JD 


27 









For our purposes, we will use the identity above at time t = T. By definition of the Bogovski operator, we have 

[ [ (p7 + $Pe)I : Vfiedxds = [ [ (p2 + 5pP)(p e - [ p e dx)dxds. (6.20) 

Jo jD Jo Jd l-k'l JD 


We can now rearrange and obtain 
r T 


[ [ Pl +1 + 5p 0 e +1 dxdt = [ [p e u e (T) ■ (p t {T) - m 0 ,s • <£ e (0)] 

Jo Jd Jd 

+ / / [2pWu e + Xdivu e I — p e u e 0 u e \ : V(p e dxds 

Jo Jd , 

rT r r rT r (6.21) 

+ / (Pi + S Pe) f Pe + fVu e V/5 e • tp e dx + / [H[div(/3 e u e )] — eV/3 e ] • p e u e dxds 

Jo Jd Jd Jo Jd 


Jo JD 
oo r T 


E 

fc=i * 


0 J D 


p e <J k (pe, PeU e ) ' <fedxd^(s). 


We proceed by estimating the p th moments on both sides of this equality. In view of Theorem IB. 51 the /? 
constraints m , and the Sobolev embedding W ¥ ' rj we obtain 


E p 


id 


(p e u e (T) - m 0 ,s)Pe(T)dx\ 


.r 

2'7 

+ \m 0 ,s\ p 2 -y 

)\pA\ 


"(W +1 

) lU 1 


2 P 

2-y 

+ m 0 ,5 2p 2 ^ ] 


if 

(U +1 ) 

r T+l 

Li x 



p i 


\ 2p 1 


Using Theorem IB.51 and (13.31) . we obtain 

r T r 


E F 

/o Jd 

< [| u e \ 2p 


2 pS7u e + Adivu e : \7ip (i dxds\ 


< E Pe [\ u e\ P ^ w i,a ) \^<Pe\ P L 2^ L 3' )i 




^H\P,\% {Li) ] 


( 6 . 22 ) 


(6.23) 


Note that (13.31) implies the embedding L& Lx P d so applying Theorem IB. 5 1 yields 


E ¥ 


< E f 


p e u e ® u e : Viptdxdsl 


o J D 
r T 


< E p 


\p e \ L p\u t \ 2 2 d | V p e \ dads 


L d 2 


L 2 P~ d 


\Pe\ 2 p\u t \ 2 2dds 

I-'X L d — 2 


<^[|Pe|^( ig) |^l3 ( ^, 2) ] 




Applying Holder yields 


E p 


/\ 

/ o Jd 


(6.24) 


(^ Pe(s)dx) JjJPz + Pl)dxds\ p < E p ‘ WIpJ^^ + \p*\l~IIz)]- ( 6 - 25 ) 


Defining r by the relation A = i + ^ — Jj and applying Theorem IB. 51 then using Holder, the embedding 

dp 

Lx f> ~ d L&, and USD, we obtain 


E f 


< E p 


0 JD 
rT 


p e u e ■ B[div(p e u e )]dxdt\ 




\p £ \ L p\u e \ _2d 7 \p e Ue\L’-dt\ 

L x 

4 p 


< E Fe 

< E S 


|p e | L g|u e | _jM_\B[div(p e u e )]\ L rdt\ 

Lx 


\Pe\ L ?J u e \ 2 ™\Pe\ L ^dt\ 


(6.26) 
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Using again the Sobolev embedding of W^ ,/3 ^ L“, we can estimate the energy correction as follows: 

r T r 


E p 

< 


/ 0 J D 


eVu e V p e ■ (p e dxdt | 

i 

For the artificial viscosity, we use Holder followed by (13.31) to obtain 


(6.27) 


E p 


ep e u e ■ V p e dxdt\ 


10 J D 


< e p/2 E Pe [|^V/5 £ |^2 (i 2 ) |/5 e |^ r(jL g ) |n e |^ (w i 




'Peru-', i 3 p 


(6.28) 




Finally, we use the BDG inequality, the summability Hypothesis 11.41 (13.31) . and the Sobolev embedding of 
L to estimate the series of stochastic integrals as follows: 


E p 


OO n'J' p 

V / / PeCTk,s{Pe,PeU e ) ■ Pedxd^(s) lP 

’0 J D 


k =1 ‘ 

00 /*T 


< E p 


^ / (/ Pe<7k,s(pe,PeU e ) ■ <f e dt\ 


k= 1 


/0 Jfl 


(6.29) 


< 


XI 

fc=l 




< E p 




l&lhl&li-* 


p/2 


p/2 




Hence, appealing to m , we may close each estimate and obtain (16.161) as claimed. 


□ 


We now proceed to a proof of the tightness. The proof is similar to the n layer tightness proof, the main 
difference being that it is less clear how to estimate the weak continuity modulus of the momentum due to limited 
uniform estimates for the energy correction term. Our strategy is to temporarily avoid this issue by proving 
simultaneously the tightness of the momentum minus the energy correction and the tightness of the momentum 
in a weaker space, then putting these together after an application of the enhanced Skorohod theorem in order 

2(3 

to deduce the tightness of the momentum in Ct ([L ] w ). Towards this end, enumerate a smooth, dense subset 

2(3 

{(/>k}k>i of [ Lx +1 } d . Define the sequence of random variables {(X e , Y’ e )} e >o via 


X r = 


(fe, Me, Peile, pi + 8p%, {/3 £ }fc> 1 


Y t = < ( p e U e (- ) , 4>k) - e / (VUeV/b, <pk)ds 


k =1 


(6.30) 


This sequence induces a measure on the space E x F endowed with its natural topology. 

„ 2/3 1,1 

E = C t (\L%) x \L1{Hq x )] w x IL?(LJ^)\ W * x [L t ;j] w x [C7 t ]°° 
F = [C t ]°° 


We may now establish the following tightness result. 

Lemma 6.5. The sequence of induced measures {P c o (X e ,Y e )^ 1 } e>0 is tight on E x F. 

Proof. We may consider each component separately, and we start by noting that the method of Lemma 15.41 
can be repeated in order to deduce the tightness of lP e o (u e , Wf) -1 , where we have denoted W e = {Pe}kLi- The 
tightness of the sequence {IP e o (p e u e , pj + 5pP ) _1 } € >o follows easily from Banach Alaogolu in view the uniform 
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bounds (EU and the integrability gains (16.161) . To treat the density sequence use the parabolic equation to 
estimate 


\dtPe\ .2 i < \PeU e \ 2/3 + C|/0 e 

Lr(Lf^) 


L™(W X ’^ fT ) 


2/3 


< , L ^ + e|p e | 




This bound can be used to show that the sequence {P £ o p e 1 } e> o assigns high probability to sets of the form 


{f&L?^ x )\\f\ , + \d t f\ _ 2 ^ < M}. 

t V *> L™(W X 


These are compact in Ct([L^.\ w ) view of Corollary IA.17I To treat the sequence {P £ o Y f 1 } £ >o, we use the 
uniform bounds (15.11) and proceed as in Lemma 15.41 to obtain the uniform estimate 


supE Pe |(/5 e u e (f) - /5 £ u £ (s),</> fc ) -e (Vh £ V/5 £ (r),<^ fe )dr| p < \<j)k\ P C i\t - s|=>. 


e>0 


(6.32) 


For any 7 < i we may choose a p large enough to ensure for each k the uniform holder estimates 


sup E F 
£>0 


1 y:\% <ic- 


We may now complete the proof by showing that the sequence of induced measures assign arbitrarily high 
probability to sets of the form 

OO 

1=1 

These are compact in view of Tychonoff and Arzela Ascoli. This completes the tightness proof. □ 

Now we are in a position to finish the proof of our compactness step. 

Proof of Provosition \6.$ . In view of Remark IA.131 E x F is a Jakubowski space. Thus, we may apply 
Theorem IA.12I to obtain a new probability space (Ha, P^), a collection of recovery maps {T £ } £ >o, and a 

limiting random variable ( Xs,Yg ) = (ps^s^pj + 5pg , Wg, Yg). Note that we may use the recovery maps to 
identify 


{peUoP? +6pP,{(p e U e (-) - f VUeVpeds,^}}^) 

Jo 

= (PeU e ,p2 +SPe,Ye)°T e -> (m S ,Pg + 5p%,Yg). 

Proceeding as in the proof of Proposition 15.31 we may obtain parts |T| and [2] along with 16.41 Our appeal to 
Theorem IA.12I immediately eives (16.51).(16.61).(16.81). and (16.91). Moreover, the uniform bounds m and the 
Banach Alagolu theorem allow us to obtain part [5] The only part requiring more justification is (16.7|l . which 
we now explain. 

Let us write Y$ = {fhg}^L 1 the limit point obtained for the sequence {Y £ } £ >o- Note that for all k and all p > 1 


lim E ®’ 5 
£->0 



r t 

p~ 

sup 

6 / (Vu e Vp € (r),</>k)dr 


te[o,T] 

Jo 



= 0 . 


Combining this with the pointwise convergence of Y e —> Yg yields the following P 5 a.s. limits 

(p e u e ,4> k ) —>• rhg in C t . 


2/3 

However, we also have the P 5 a.s. convergence p e u e —> pu in )]„,*, which identifies a candidate limit 

at the outset. Arguing similarly as in the corresponding proof of Proposition 15.31 we may obtain (16.71) . □ 
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6.2 Preliminary Identification Step 

We now define the filtration {J r t 5 }, 5 >o using the restriction operators with respect to the spaces dictated by 
Pronostion 16.31 

Tl = <j(r t ps , r t us, r t ( pgug ), r t Ws ). (6.33) 

Lemma 6.6. The pair ( ps,ug ) satisfies the continuity equation. \6. c 2\ of Defjnition \6.‘A 

Proof. In view of part [2] of Proposition 16.31 the pair (p £ ,w £ ) satisfy the same parabolic equation on the new 
probability space (Qg, Tg, Pa). Hence, in view of (16.511 and (16.71) . we may pass e —> 0 in the weak form of the 
parabolic equation for ( p e ,u e ) to obtain the desired result Pa almost surely. Indeed, the only additional detail 
is to argue that the viscous regularization tends to zero. Recall the energy identity used in the n —> oo density 
upgrade. In view of the representation of (p e ,u e ) in terms of the recovery maps, this identity holds on the new 
probability space, so estimating the moments on both sides(neglecting the initial data, which are controlled) 
yields 


Iv^Vp, ' 2p 


e 'L 2 t (Ll) 


< E P5 


| p £ divi 




< 


2 p |~ IP 


< 




12 p 


Hence, we may use the uniform bounds obtained in (15.11) and estimate for <fi £ Cf°(D) 


\S7 p t \\S7 cj)\dxds\ p 


/ o J D 


< e%\S7(j)\ L iWF s \\y/eVp, 


elj dUj 


0 . 


□ 


Lemma 6.7. The pair ( ps,ug ) satisfies the momentum equation (16.31) from Definition \6.Sl with a pressure law 

Ps + S Pg ■ 

Proof. Let us proceed with the method in Lemma ISUl bv defining for each <j> £ C°°(D) the {Jy 5 }^1 0 continuous, 
adapted stochastic process {Mf (<p)}J =0 via 


Mf ((/))= / pgug(t) ■ <f> dx — / m 5 0 -(j)dx — 


D 


D 


n [p5Ug ( 


) us — 2fiVus\ : Vfidxds 


+ f f \{p$ + 5pg — Adivw^)/] : \7<j)dxds. 
JO J D 

Note that by appealing again to (15.Ill , we have 


E Pi 


eVii £ Vp £ • (j)dxds\ 


0 J D 


< \4\l° 


12 p 

el L?(Wf 2 ) 


l\/eVp £ | i ^ (i 2) 


—> 0 . 


Combining this observation with Proposition l6.3l vields all the necessary ingredients to proceed with the method 
in Lemma 15.61 and identify 


OO p-fc p 

m) = E / / Ps a k{ps,Psug) ■ <j>dxdfig{s). 

k=i Jo Jd 


□ 
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6.3 Strong Convergence of the Density 


Now to proceed to the proof of the strong convergence of the density. The first step is the following weak 
continuity result: 


Lemma 6.8. Let K CC D be arbitrary, then the weak continuity of the effective viscous pressure holds on 
average, that is: 


lim E*" 5 [ J J (( 2 /r + A) divu e — pf — 5 pf) p e dxdt\ 

T r _ _ 

((2/i + A) di vug — p} — 5pg)psdxdt}. 


o Jk 


Proof. Recall that A = VA" 1 , where the inverse laplacian is understood to be well defined on compactly 
supported distributions in R d . Let 77 be a bump function supported in D. Define the following two random test 
functions: p t = rjA [r]Pt] and ps = rjA [rips]. Using the parabolic equation for p e driven by u e and the transport 
equation for ps driven by us, we may check that 


d t pe = rjA o div(r](eVp e - p e u e )) + r/A [V 77 • ( p e u e - eVp e )] 
dtps = —r)A o div (p psus) + pA [V 77 • psus] 


Using the momentum equation for ( p e ,u e ) and ( ps,us ) we may use the Ito product rule twice(see the remarks 
in Proposition 16.41 regarding justification) to find the evolution of p e u e ■ p e and psus ■ ps- The first application 
yields the P 5 a.s. equality 


/ p e u t (T) ■ p e (T)dx = / 77io,,5 • Pe{0)dx + / / [p e u e ■ d t p e + [peU e ® u e — 2p\7u e ] : \7p e ]dxdt 

J d Jd ’ Jo Jd 

+ f [ ( — Adi V U e + pf + 5pP)I] : Wpe] — e'VUe'V p e ■ p e dxds 

Jo Jd 

00 „fp „ 

+ Y] / Pe<?k,s(Pe,PeU e ) ■ Pedxdf3s(s). 

k= i J o Jd 


The second application yields the P 5 a.s. equality 

/ psus(T) • (ps(T)dx = / m 0 ,s • <ps(0)dx + / / [psus • d t <ps + [psus 0 us ~ 2pS7us\dxdt 

Jd Jd Jo Jd 


+ / / [(—Adi vus + p] + 5pg)I] :V<j>8\dxdt 

Jo Jd 

OO njf n 

+ V / / pscr k ,s{ps,psus) ■ psdxd/3g(s). 

k =1 Jo Jd 


Note that 


[p e u e <g) u e — 2 pVu e + (—AdivM e + pf + 5p^)I] : \7p e 
= [p e u t ® u e — 2 pX/u e + (—A div u e + pf + 5p @)/] : V?; (g> A [pp e ] 

+ pp e u e <X> u e : S7A[pp e ] — 2ppVu e : VA[pp t ] + p(pf + 6p f — Adivu e )p e . 


Moreover, integrating by parts twice(justifying on a smooth approximation) reveals 

f f -2ppVu e : S7A[pp t ] = f [ 

Jo Jd Jo Jd 


2pp 2 div u e p e dxds 


Also note that <^ € (0) = ^$(0). Taking expectation(so that the stochastic integrals vanish) of both Ito product 
rules above yields two fundamental identities: 


f f r] 2 [(2fi + A) divit e — pi — 5p^]p e dxds 
Jo J D 


= / U 


'0 J D 

, tA 


(6.34) 


r c,e 


T C,e 


T C,e 


f f - 


7 ? 2 [(2 p + A) divua — fig — 8p$]pdxds 


= r + /; 


c 


I? 


I3 


l 2 ■ 


(6.35) 


Our labeling convention should be interpreted as follows. The terms d) 4 ' 6 ,/^ 4 ’ 6 are “artificial ” and will tend to 
zero as e —> 0, If"’ 6 ,I 2 e , are lower order cutoff’ terms arising due to the localization of the estimate, and 
I 1 ’ , / 2 ’ are the principal terms arising irregardless of the boundary conditions. More precisely, the contribution 
at the e layer yields 


r = : 


If' e =: 


i£’ e = —i 


/f - 6 = E ^ 4 


Ud 


rimo,s ■ A[ppo,s\dx 


o Jd 

rT 


erjp e u e ■ A[d\v{rfJ p t ) — V?y • V p e ]dxds 


eVUeVpe • p e dxds 


io Jd 


n 

/o Jd 


+ : 


=: 


io Jd 

rT 


[p e u e 0u e — 2pVu e ] : Vr? <g) A [pp e ] dxds 
(—A div u e + pi + 5pP)I] : X7p 0 A [pp e ] dxds 
p t u t • M[V?7 • p e u e ]dxds 


lo JD 


lg’ e = E ®’ 5 


A p ’ e =: 


J 2 P ’ £ = -i 


to Jd 
r T 

/ 0 JD 


u e ■ [A{jpp t )S/p — Vppp e ]dxds 

/ p [p e u e 0 u e : X7A[pp e ] — p e u e ■ A o div(^p e u e )] dxds 
Jo Jd 

p e u e (T) ■ <p e (T)dx 


id 


(6.36) 

(6.37) 

(6.38) 

(6.39) 

(6.40) 

(6.41) 

(6.42) 

(6.43) 

(6.44) 
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In the limit as e —> 0 we expect to obtain the following contribution 


if = : 


+ ] 


I 2 = 


n 

I o Jd 


/0 JD 
r T 


Jf = E p ' 5 


/f = E Ps 


Lf = -E P4 


/o Jd 

r T r 


[psus — 2 pWus] : Vr/ <g> A [77/35] dxds 
(—Adivu^ + / 3 ^ + fips)I] '■ V77 <S> -4 [77/05] dxds 
P5U5 • Al[V?7 • pswajdxds 
m, 5 • [Al( 7 ?/ 3 i)V ?7 - \7r)r]ps\dxds 


) 0 JD 
r T r 


I / 77 [/95M5 <S) Mi : V.4 [77/05] - psus ■ A O div(rjpsus)] dxds 

0 JD 

p S u S (T) ■ (ps(T)dx 

L JD 


(6.45) 

(6.46) 

(6.47) 

(6.48) 

(6.49) 

(6.50) 


First note that by the uniform bounds in (15.11) combined with an interpolation argument we obtain the estimate 


it'* + J 2' e ;$ [| v^Vp £ 


| PeU e 


\ U ^Ll(Wl' 2 )\ 


0 . 


(6.51) 


To treat the remaining integrals, note that by the uniform bounds (15.11) and the Vitali convergence theorem, it 
suffices to establish the relevant Pi convergence, so the analysis essentially reduces to the same arguments as 
in the deterministic framework(by design). We recall them here for the convienience of the reader. 


Starting with the first cutoff term, note that for all q > 1, A : L% —> L%. is compact. Hence we may com¬ 
bine (16.51) with appendix result Theorem IB.51 to obtain A[r/p e ] A[rjps] strongly in L\ x . In view of a similar 

fid 

argument in Lemma f5.61 we have p e u e ® u e —> psus ® us in L 2 {Lx'' d ~ 1)+d ). Combining these two observations, 
(16.61) and (I 6 . 8 |) we end up with a product of a weakly converging sequence and a strongly converging sequence 
and conclude I^’ e —> Xp. A similar argument also yields d p,e J p . 

2/3 

Next note that A : L£ +1 —► L x is compact for i > \ We may now use (16.71) together with 

Theorem Ib31 to conclude that A[Vr/ ■ p e u e ] —> A\S7p ■ psus] strongly in L™(L X ) for all m > 1 and r satsifying 
the relation above. We may use again (16.71) to obtain another weak times strong where the exponents match 
up appropriately since /3 > d, allowing us to conclude I — y . To treat the final cutoff term, simply argue 
as in the passage to the limit in the flux term of the continuity equation and obtain / 3 ’ —> Hr. 

The treatment of the principle term d ^ 3,1 is the nontrival part, but we have built most of the work into an 
appendix result based on the Div Curl lemma. Working componentwise we may write 


l[’ e = y E Ps 


i,j =1 


’ D 


ul(rip e dij A 1 (r/p e u J e ) — r/p e u J e dijA 1 (pp t ))dxds 


(6.52) 


In view of (EH and (1CT7D . we may appeal to Lemma IB.41 with p = /? and q = , making use of the compact 

2/3 , 

embedding ^ W x ’ for p > |d in order to conclude that P 5 a.s. 

(rjpedijA-'trjpeul) - pp t u{d i jA- 1 (rip t )) (ppsd lJ A~ 1 (r]psu 1 s ) - T}psu 3 g d i jA~ 1 {r}psj) 
strongly in L 2 (W~ 1,2 ). Appealing once more to (15.11) and Vitali we find df ’ 6 —> d-f. □ 

We now proceed to a proof of the strong convergence of the density. 
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Lemma 6.9. The sequence of densities {p e } e >o converges strongly to ps in the sense that for all p > 1 and 
r < /3 + 1 


lim |p e - p s \ 


e->0 " ' 'Lp(p s -,L £ x ) 


= 0 . 


(6.53) 


Proof. . We will begin by establishing pslog(ps) = pslog(ps) almost everywhere in Qs x [0,T] x D. Using 
the Di Perna Lions commutator lemmas and the Hardy inequality, we may renormalize the parabolic equation 
for the p e and the transport equation for ps separately (and note they have the same initial data), using the 
renormalization /3(p) = plog(p). As a result, we find that for all smooth if with if{T) = 0 the a.s. inequality 
holds 

/ 1ft [ps log ps - Pe log Pejdxds < if [ps div Us - Pe div Ue ] dxds. 

Jo Jd Jo Jd 

Let us further decompose the RHS into a portion over a compact set K CC D where the weak continuity result 
Lemma 16.81 may be applied, and a remainder R e K {if) 


/ / ift[ps log PS - Pe log Pe]dxds < / / if [ps div Us - Pe div Ue] dxds 

Jo Jd Jo Jk 

+ / / if [ps div us — Pe div u e ] dxds 

Jo J D\K 


10 JD\K 

r T r _ — 

ifps(divus — Ps — Sps) - tfpe(divu e - pj - 5pP)dxds 


> 0 J K 
rT 


10 JK 


[ ifPsiPs + 6 Ps)~ i’PeiPl + &Pe)dxds + R^lf). 

Jk 


We now take expectations with respect to on both sides of the inequality above and send e —> 0. By Lemma 
(16.81) . the first term tends to zero. Moreover, we may apply (16.131) to conclude the remainders R e K (if) converge 
to Rx^if) defined below. Hence, we conclude 


fj 


iftE Ps [p s log ps - ps log p s \ dxds 

[ i>Ps(Ps + s Ps)- i’peipl + 5p^)dxds 
Jk 

Rxi^f) = E 1 ” 5 I j if [ps div us ~ Ps div dxds 


< lim inf E Ps 
€->0 


+ (VO¬ 


ID J D\K 


(6.54) 

(6.55) 


Applying Minty’s trick to the convex function fj{p) = p 1 + 5pP , the first term on the RHS of (16.541) is non¬ 
positive. Hence, for any lebesgue point s € [0,T] of the function s E p,s log pa — pa log pa] (s), we may 
choose a sequence of test functions that approximate l[o, s ](f), so that their time derivatives approximate the 
negative of a dirac mass centered at the point s. Hence, the following inequality holds almost everywhere in 
time 

77 * 

Jo J D\K 

We can make the RHS of (16.561) arbitrarily small by choosing K close enough to D. Hence, the quantity on 
the LHS of (16.561) starts from zero, is always positive(in view of the convexity of plog(p)), and almost never 
increases. Thus, the LHS vanishes almost everywhere in time and we conclude p^log^) = ps\og(ps) almost 
everywhere in Oj x [ 0 ,T] x d. 

Combining this with the weak convergence (16.141) . we may apply Lemma rB.3l to conclude the ilj x D x [0, T] 
almost everywhere convergence of p e ps away from the vacuum regions of ps. Moreover, since {p e } e >o is a 
nonnegative sequence, weak convergence implies almost every where convergence on the vacuum regions of the 
limit. Finally, we may use the uniform bounds (15.11) in order to deduce the desired claim. □ 


[ps div us — ps div dxdr 


(6.56) 


E P5 


{ps log ps - ps log ps)dx 


L JD 


(s)< 
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6.4 Conclusion of the Proof 

Proof of Theorem \6.1\ For each S > 0 we apply Proposition 16.31 to construct a candidate S layer approximation 
( pg , ug). Combining the preliminary limit passage, Lemmas 16.61 and 16.71 together with the strong convergence of 

the density, Lemma 16.91 we are able to identify pf + Spg = pj + Spg and complete the identification procedure. 
The uniform bounds can be argued as in the proof of Theorem 15.11 □ 

7 Proof of the Main Result: 5 —* 0 
7.1 S —> 0 Compactness Step 

Following Feiresil [TO] , we introduce a collection of approximations from below of the renormalizations /3(p) = p 
and /3(p) = plogp, respectively. These will be useful in our proof of the strong convergence of the density later 
on in the section. Namely, let T(p) = pi[o,i](/°) + 0(p) 1 [ 1 , 3 ] (p) + 2 l[ 3 j0 o)(p); where <f> is chosen to ensure global 
smoothness and concavity of T(p). Use T to define the approximations by letting for each fc, T k (p) = kT(^) 

and L k (p) = p ff 1 k z [' 1 ■ Note that by design pL' k (p) — L k {p) = T k (p). We now proceed to our final application 
of the Skorohod theorem. 

Proposition 7.1. There exists a probability space (U, T, P) along with a sequence of “recovery” maps {Tjsjax) 
and limit points 

(p, u, pT {/3 k }%Li, VP U , plogP, {Tc(p)}^ 1 ) 

({W}r=i,m(p)divn}r=i, 

f s -.(Sl,F,V)->(Cl s ,P s ,h) 

such that the following hold 

1. The measure P 5 may be recovered by pushing forward P under Tg. 


2. The new sequence {(p, 5 , Ma)}, 5 >o defined by (pg,ug) = ( ps,ug ) o Tg constitutes a S layer approximation 
relative to the stochastic basis (Cl,J 7 ,¥,{J r g}t>o,Ws), where 

Wg := Wg o Tg Tg := a ( r t ps , r t u S: r t (pgug), r t Wg) 


3. The following uniform bounds hold for all p > 1 


sup Hr 
5>0 


+ + !«»«!" (i2l 


4- Let q £ ( 1 , 7 ), then the following convergences hold P a.s. 

PS -t P 
ug —^ u 

pgug pu 

P's ^T 7 

Wg^W 

Tk(ps) —> T k (p ) 

Lk(ps) —> L k (p) 

(. PsTk(pg ) - T k {pg)) div ug -> ( pT' k (p) - T k (p)) div u 
ps log(pa) plog(p) 


+ \ u s\ 2 A 


(K 


in 

in 

in 

in 

in 

in 

in 

in 

in 


) < 00 

(7.1) 

c t {[L2\ w ) 

(7.2) 

L t( H 0,x) 

(7.3) 

C t {[Lp] w ) 

(7.4) 


(7.5) 

[C t ]°° 

(7.6) 

4 'yd 

C t {[L^] w ) 

(7.7) 

c t mu 

(7.8) 

L?{Ll) 

(7.9) 

c t mu 

(7.10) 
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5. The following additional convergences hold: 


yfFws —t \fP u 

m Ll^L?{Ll)) 

(7.11) 

PS -> p 

m L^(n;L?(LZ)) 

(7.12) 

PS log (pa) -t plog(p) 

in If w .(Sl-,L?{LD) 

(7.13) 

Tfe(pa) div us -»■ T k (p) div u 

in L p w {Ll x ) 

(7.14) 



(7.15) 


In order to prove the tightness of the pressure, we will need the following integrability gains. 
Proposition 7.2. The following integrability gains hold for arbitrary p > 1 and n < min(^ — 1, ^): 

r T 


supl 

<5>o 


/ [ Pl +K + 5p 0 s +K dxdt\ 

Jo J D 


< oo. 


(7.16) 


Proof. The proof of this proposition follows along the same lines as the proof of the integrability gains at the e 
layer. We will give a limited amount of details and mainly indicate how the constraint on k arises. Introduce 
the following random test function 

<Ps =B[Ps ~ j^j- f D Ps dx \- ( 7 - 17 ) 

By the Di Perna Lions commutator lemmas, we may check that Pa almost surely, in the sense of distributions 

1 


Ps div(wa) - 


\D\ 


Pi div(us)dx 


= 0 . 


9t{Ps - ]^| Ps) + di v{ffiu s ) + (k - 1) 

Applying the Bogovski operator on both sides yields 

d t <ps = -B[div(paua)] - (k - 1 ) 1 ? p's div(ua) ~ J Ps di v{u s )dx 

We may now use the Ito product rule together with the momentum equation to find the evolution of psus ■ fis, 

[ [ Ps +1 + SPs +1 dxdt = f [p s u s {T) ■ <p&(T) - m 0 ,s ■ <P<s(0)] 

Jo j D Jd 

+ / / [2pVus + \ div u$I — psus 0 fi*] : V<ps 

Jo Jd 

+ [ [ (Ps+SPs) i TFu dx + f f psus ■ B[dW(tfus)]dxds 

Jo Jd Jd \ u \ Jo Jd 

+ (k ~ 1) / / P 5 U 5 * B p's di v(u s ) - f Ps di v(u 6 )dx 

Jo Jd L Jd 


dxds 


E 

k= 1 


0 JD 


pscrk,s{ps,psus) ■ psdxdps(s). 


We estimate the new term from the renormalization as follows 
E*' [|pm-B(p?divu«)|£ 






1/3 


M 


3m 




1/3 


E Ps 


\B(ftdWu 5 )\% {L 


1/3 


where 1 = 1 + ^ — i. Define the exponent p via ^ | — L. The Sobolev embedding Wf’ p ^ L r x gives 


\B(ftdivu s )\% (L 


< 


\PS\^(LI) 


t \ X / 
- 11/6 


1/3 


< 


< 


I Ps div««liT (IiS) 

1/6 


1/3 


divusl^L*) 
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where ^ ^ ^ by Holder, which is valid in view of 7 > tj. To control this final term we require 

Kq < 7 which leads to the condition k < -J — 1. Similar estimates yield uniform control of the terms 




\J3[div(p%us)] ■ p s us | 


and 


\p 5 u 5 ® us : '\7B{pf)\ 1 ff 1 . The stochastic integrals and the term 


E Ps 


\psus(T) ■ MT)\Ti 


are both controllable under the condition k < ^ — 1. Finally, the dissipative term and mean value correction 
can be estimated provided that k < 

Define the sequence of random variables {(Xa, Ya)}( 5 >o via 
X s = (ps,us, psus, pg, ps log ps, 0s}kLi) 

Ys = ({T.ips^ALkimkLidiPsTUM -T k (p s ))dWus}r=i). 

These random variables induce a measure on the space E x F where 


□ 


(7.18) 


E = C t (\Ll \ w ) x L 2 t (Hl x ) x C t ([LS +1 U ) x C t ([L%] w ) x [l£*] w x [C t ] 


F = 


d 

■ 2 ^ 


C t {[L?- d ] w ) x C t {[Ll\ w ) x Li iX 


(7.19) 


The space is understood to be endowed with its natural product topology. 

Lemma 7.3. The sequence of induced measures Pa o (Xa,Va ) —1 are tight on E x F. 

Proof. We proceed componentwise. Note that the methods in Lemma 16.51 can be used to deduce the tightness 
of {Pa o (^ps,us, pj, Ws'j }a>o- To treat the momentum sequence, argue as in the other layers and deduce for 
all cj) S C x the following estimate holds 


swpE Ps \(p s u s {t) - psus(s),(/))\ p < \<j)\ v cl \t- s|S. 


<5>0 


(7.20) 


In view of Lemma TA. 181 the tightness of {Pa 0 (/5atta) _ 1 }a>o follows. Next we will treat the induced measures 
Pa 0 Ys- By the Di Perna Lions commutator lemmas, the following equality holds(in the analytic sense of 
distributions) Pa a.s. 

d t L k (p s ) = - dw(L k (ps)us) -T k (ps)dwus- 
Noting that L k (p) < kp and T k (p) < k we obtain the estimate 

\d t L k (p s )\ 


-1, A2 

"7 + 2 * 


^ |ifc(P5)lL“(L2)I^Lj(L2*) + \Tk(ps)\L™ x | divuaL^ 


L?(W X ) 

^ (&+ \ps\l™(lZ))\us\l-}(hI x ) 

A similar estimate holds for T k (ps). Let {M k }’^L 1 be a sequence of real numbers and note that 


_ 4 'yd 

G C t ([Lf^] w ) x C t ([Ll\ w ) x L\ x \ \d t f\ 


k =1 


L 2 t (W x 


I/I 




4 "yd 


+ \d t g\ 


L 2 t (W a 


-i. 

’ 7+2* 


\g\L^(Lg) + 1 < M k 2 k } 


is a compact set in 


4-7 d 

• ^7 


C t ([Lj^U x C7t([ig] w ) x 


by Corollary 1A. 171 , Banach-Alagolu, and Tychonoff. 


Hence, we may appeal to (ED and Chebyshev to show that sets of this form have Pa probability arbitrarily close 
to one, uniformly in S. Finally, one can treat {Pa 0 {ps log/5a) _ 1 }a>o a similar way, using the renormalized 
form of the continuity equation and the integrability gains to control the term ps div us- □ 


Proof of Proposition \ 7.1\ : The proof follows along the lines of the e layer. 


□ 
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7.2 (5—^0 Preliminary Limit Passage 

Using the restriction operator according to the spaces described in Proposition 17.11 we define the filtration 
{J r t}T= o via 

Ft = a(r t p,r t u,rt(pu),rtW) (7-21) 


Lemma 7.4. The pair ( p,u ) satisfies the continuity equation, part \l. 7| of 1 1 . 51 

Proof. In view of the strong convergence of the initial density laid forth in Hypothesis 11.21 we may proceed as 
in Lemma 16.61 □ 

Lemma 7.5. For all <j> £ Cf°(D), the process defined by 

M t ((f>) = / pu(t) ■ fidx — 

J D 

is a continuous, {.F 4 }^Lq martingale satisfying for all p > 1 


/ mo • </> — / / [pu ®u — 2/.tVu — A div ul\ : X7cf> + p 1 div fidxds. 

J D Jo Jd 


E 


sup \M t (4>)\ p 
te[o,T] 


< oo 


(7.22) 


Proof. Introduce the continuous process {Mf (</>)} f =0 defined by 


(fi) = / p s u s {t) ■ fdx - / m 0 ,s •</> — / [psus ® ug - 2p\7u s - A div :V(f> + p] div 

Id Id do Jd 


fidxds 


We may use Proposition 17.II along with Hypothesis 11.21 to establish for all t £ [0,T] the convergence —y 

almost surely with respect to P. Indeed, the only additional steps(other than what was required at the 
e layer) are noting 


E 


/ 0 JD 


Spsdx l J 


,8p 

< 6/3 + *E 


f j pd +K dxds\ 

Jo Jd 


0. 


and also that the strong convergence of the intial data holds by Hypothesis 11.21 The estimate above also leads 
to the uniform bounds 


supE 

<5>o 


sup I M t s (fi)\ p 
te[o,T] 


< oo. 


(7.23) 


This information is enough in order to use our usual procedure and verify that {Mt{<j>)}t> o is a continuous 
{.F 4 }^ martingale. Additionally, we may check the convergence in (fi; L°°[0, T)), so the estimate 17.221 
follows. □ 


7.3 Strong Convergence of the Density 

We now want to work towards establishing the weak continuity of the effective viscous pressure. At the e layer, 
we chose a test function and begin by applying the Ito formula to find the evolution of psu$ ■ ip$. Since we have 
not identified our martingale as a stochastic integral, this is slightly less straightforward. Instead, we check only 
that our desired identity holds in P expectation. Let us now define the following two random test functions: 


<Ps(t,x,w) = i]A[pT k {p 5 )} 


p(t,x,uj) = r]A yr)Tk(p) 

By the Di Perna Lions commutator lemmas, we may verify the following identity: 

dtps = -pA o di v{r]T k (ps)us) + rjA [V 77 • T k {ps)us] - 77 A[r)(psT' k (ps) - T k (ps)) divwa] 
Sending <5 —>• 0 and using the Skorohod step we find 


d t ip = -r)Ao &vvfqT k (p)u) + 77^4 V?j-T k (p)us - gA[i](pT^(p) - T k {p)) div it] 
We now establish the following averaged Ito product rule. 


(7.24) 

(7.25) 
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Lemma 7.6. (Two Averaged Ito Product Rules) Define 4>s und <f> as above, then the following two averaged Ito 
product rules hold P a.s. for all times t € [0,T] 


E p 


-E* 


E p 


-E* 


pgugft) ■ <psff)dx = / m 0j 5 • ips(0)dx + E p / / \pgug ®ug — 2 p\7ug] : Vp$\dxds 

JD ' JO JD 

n {(—\div us + p) + 6pg)I] : V(ps\dxds+ E r f f [psus ■ d t <ps]dxds. 

) Jo Jd 

puff) ■ ip{t)dx = / mo • tp(0)dx + E p / / [pu ® u — 2 p\7u] : \7<pdxds 

d J Jd Jo Jd 

n [(—A div u + p 7 )/] : V<p + pu ■ dtipdxds. 


Proof. The first identity can be proved in the same way as at the e layer. To proceed to the second, denote by 
the standard mollifier(localized at scale k) and f KtX the standard mollifier centered at the point x. Extending 
(p, u, Tk(p ), p 7 , pP) by zero outside of D, we may define the quantities g K , p K via 


g K (x,s)= / \pu® u(s) - 2p\7u(s) + (-X div u(s) + p^(s))I] \ \7f K (x - y)dy 

J R d 

<Pk(x, t) = (<p(t) * ^k) (x) (pu) K (x, t ) = {puff) * f K ) (x) 


By definition of the process {M t {f K ^ x )}t >o we obtain for all k > 0 and x £ D, the following equality holds P a.s. 


(pu) K {x,t) = m 0}K ,{x) + / g K {x, s) + M t {£ K , x ) 

Jo 

By Lemma [7751 the process {M t {f K}X )}’(' =0 is a martingale satisfying enough bounds to give a meaning to the 
stochastic integral below. Applying the classical Ito product rule for continuous one dimensional martingales, 
we obtain for each x £ D, the following equality holds P a.s. 

(pu) K (x, t) ■ <p K {x, t) = m 0iK ,{x) ■ <p K (0) + / [<p K {x, s) ■ g K {x , s) + d t ip K ■ ( pu) K {x , s)] ds 

Jo 

T / Tk{x, s)dM s {f K ^ x ). 

Jo 

Note the estimate 


E 

[ c/) K (x,s) 2 dM(^ x ) s 

< E 

sup \<t>k{x, s) 4 


Jo 


i£[0,T] 


1/2 

E \(M{Z k , x ))t\ 1/2 ■ 


By the definition of quadratic variation(or the Doob Meyer Decomposition for continuous martingales) and the 
uniform bounds on the fourth moments of Mt(£ KjX ), the second moment of the quadratic variation is controlled. 
The other term is estimated using the expression for dt4>k,x implied by the equation above. Hence, the stochastic 
integral above is a martingale (rather than just a local martingale), and hence has mean zero. Taking expectation 
and integrating over D yields 


E p 
+ E 1 


{pu) K {x,t) ■ ip K {x,t)dx 


= / mo,K.{x) ■ ip K {0)dx 
Jd 

/ / [tp K (x,s) ■ g K (x,s) + d t <p K ■ (pu) K (x,s)]dx 


L^o Jd J 

Letting k —> 0 and appealing to standard properties of mollihers, we obtain the result. □ 

We may now use our averaged Ito product rule, together with the Skorohod step in order to obtain another 
weak continuity result. 
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Lemma 7.7. Let K CC D be arbitrary, then the following averaged version of the weak continuity of the 
effective viscous pressure holds: 


lim E p [ j [ {p] + Sp f s - (2/z + A) di v us)T k (ps)dxdt\ 

S ^° Jo Jk 


(7.26) 


= E F [/ / (p~t — (2p + X) div u)T k (p)dxdt\ 

Jo Jk 


Proof. In view of Lemma l7.61 we have again two fundamental identities the drive the result: 


E p 


f [ p 2 [(2p + A)divu 5 - p] - 6p%]T k (p s )dxds 
Jo J D 


10 JD 

= I 0 ’ 5 + if ’ 5 + if ’ 5 + if ’ 5 + I R ’ S + I P ’ 5 

r T 


E F 


[ I ? 7 2 [(2/x + A) divu — p' 1 — 6pP]Tk(p)dxds 

Jo Jd 


= 1° + if + if + if + I R + I P 


Our labeling convention follows the same logic as the decomposition in Lemma |6.81 There is only one term of 
a new character, I R ' S arising from the renormalization of the continuity equation by T k (p). 


I°’ s = / pm 0 ,s ■ A[pp 0 ,s]dx 
Jd 


if ’ 5 = E f 


if ’ 5 = E f 


i 3 c ’ 5 =: 


I R ,s = E P 


J p ,s = E P 


/ / [psus 0 us - 2pVus + (-Xdivus + p])I] : Vp 0 A [pT k {ps)\ dxds 

Jo J D 

psus ■ A[S7p ■ T k {ps)us]dxds 
u s ■ [A{pps)Vp - S/ppps] 
ppsus ■ A [p(p s Tf(ps) - T k {p s )) div 14 ] dxds 


lo JD 

r 

0 J D 

fl 

Jo Jd 

r T r 


Jo JD 

f 

Jo Jd 


p[psus 0 us : S7A[pT k (p s )\ - psus ■ A[div(pT k (p s )us)}]dxds 


I°= pm 0 ■ A[pp 0 ]dx 
Jd 


if = E p 


if = E p 


if = E p ' 5 


I R = E f 


i p = e f 


/o Jd 

f 

Jo JD 

rT r 

Jo Jd 

rT r 


[pu 0 u + 2/.tVu + (A div u — p 1 )!] : Vp 0 A pT k (p) 


dxds 


pu ■ A[Vp ■ T k (p)u]dxds 


u ■ [A(pp)Vp — Vppp\ 


10 Jd 

f!' 


ppu-A p{pTf.(p) - T k {p)) div u 


dxds 


pu®u : VA[pT k (p)\ - pu ■ A[div(pT k (p)u)\ 


dxds 


(7.27) 

(7.28) 

(7.29) 

(7.30) 

(7.31) 

(7.32) 

(7.33) 

(7.34) 

(7.35) 

(7.36) 

(7.37) 

(7.38) 

(7.39) 
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Note that too, a —> too strongly in L ^ TT , hence it suffices to note that Wf n c —> in order to conclude that 

I 0,s —> 1°. The remainder of the analysis is devoted to working term by term and showing the convergence of 
each integral, similarly to the analysis at the previous layer. This is accomplished by combining the information 
in Proposition 17.11 with appendix Lemmas IB. II and IB.41 □ 

Let us proceed by applying weak continuity result to establish the strong convergence of the density. 

Lemma 7.8. The sequence of densities {ps}s>o converges strongly to p in the sense that 


$m\Ps ~ p\ LP (n-,L*-) =°- 


(7.40) 


Proof. Our strategy is the same as in the the previous layer. The hypothesis m ensures that 7 is large 
enough that the integrability gains imply p £ L p (fl;L^ x ) for some p > 2. Hence, weak solutions to the 
continuity equation are still renormalized solutions. That is, it remains valid to apply the commutator lemmas 
and the Hardy inequality and obtain for each smooth ip with ip(T) = 0 the following equality holds P a.s. 




ipt[L k (p) - L k (p s )\dxds 



divit — T k (ps) div ps\ dxds. 


In fact, the equality above should also contain the initial data, but since they are converging strongly in view 
of Hypothesis 1 1. 21 we simply take the data to be zero. Decompose the RHS into a portion over K CC D and a 
portion near the boundary 


/7 


>0 J K 
pT r 


il)t[Lk(p) ~ Lk(ps)\dxds = 

Jo J K 

[ ^ T k (p)(divu-~fp) -T k (ps)(divu6 - pj) 

dxds 


ip [Tk(p) div u-Tk(ps) div ps) dxds + R 5 K (ip) 
dxds 


10 Jk 


ip T k (p) p" 1 — T k (ps)p1 


Taking expectation with respect to P on both sides and using the weak continuity result Lemma 17.71 we obtain 


E 

/ / ip t [L k (p) - L k (ps)\dxds 

< liminf E 

/ T / 

Tk(p) p 1 -T k (p s )p] 

dxds 


Jo JD 

(5—>0 

Jo Jk 




- E 


la Jk 


ip[T k (p) - T k (p)\ div u dxds 


+ 


/ 0 JD\K 


ip T k (p) div u — T k (p) div u 


< | divu| I/ 2(Q. i 2( i 2))|T fc (p) — T k (p) |i2(n ;i 2 j 


0 Jd\k 


ip T k (p) div u — T k (p) div u 


dxds 

dxds. 


Note that we used the montononicity of the pressure. The uniform bounds on the density in L p (f2;L^ a .) 
imply that the first term tends to zero as k —> 00 (by interpolation and a straightforward lower semicontinuity 
argument). Hence, using a sequence of ip approximating indicators(as in the e layer) and making the remainder 
term arbitrarily small, we find by taking k — > 00 that for all t £ [0, T] 


E[ / plog(p)(t) — p\og(p)(t)dx\ =0. 
JD 


Arguing as we did in the previous layer and using the uniform bounds, we conclude. 


□ 


The following lemma is due to Lions ED, and will be used to deduce the strong convergence of the momentum. 


Lemma 7.9. Let {pn}ffL\ and be deterministic sequences satisfying the following 
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1. For all p <7, p n -> p in C t {L p ) n C t ([L].] w ) 

2. u n —u in L 2 {Wf’ 2 ) 

3. pnUn^-pu in Ct([LS +1 ] w ) 
then p n u n —► pu in L 2 {L X ). 

Proof. See [2T] page 34-35 for most of the ideas of the proof. The rest are left to the reader. □ 


7.4 Conclusion of the Proof of Theorem 11.61 

Lemma 7.10. The pair ( p,u ) satisfies the momentum equation (11.811 of Definition \1.5[ 

Proof. We proceed with our usual strategy based on Lemma TA. 151 Namely, for each f> £ C£°(D) we introduce 
the continuous, {.F 4 }^o adapted process {M t }^ 0 defined by: 


Mt{<jf) = f pu[t)-<f)dx— f mQ-cfdx— f f [pu ® u — 2pS7u — Adivit/] : V</> + p 1 div fidxds . (7.41) 

Jd Jd Jo Jd 

In view of Lemma m and Lemma [7751 . {M t (4>)}J =0 is a martingale. To complete the usual strategy, note that 


E 1 ^ 


/ 0 JD 
rT 


< e f 


■ e ¥ 


( J \ps°k{ps * Vs,{psus) * VS,x) - pcr k {p, pu,x)\dx) 2 ds 
I ps - p\Wk{ps * Vs, {psus ) * vs,x)dx\ 2 ds 
(J p\<J k {ps * Vs, (psus) * vs,x) - a k (p, pu, x)\dx) 2 ds 


' 0 JD 

rT 


IQ JD 


< IrrJ 2 


O’ki 2■ 


-2X. 

(L“ ) 


E 1 ^ 


I PS ~ P \ 2 _2"Y 


-E f 




E 1 




ry 2(1 _ ‘2 

I <Jk{ps * Vs, {psus) * vs,x) - <j k {p,pu,x)\^dx) 7 ds) z 


I 7 \ 2(1 — —■) . ^2 

\Uk[PS * vs, {PSUS) * VS, X) - Ok\p,pu, 'X) — ~ 

>0 Jd 

The hrst term tends to zero by Lemma T7.81 To control the second term, note that 

rT 


1/2 


r 


E- 


< E 1 ^ 


(/ (/ Wk{ps*Vs,{psus)*V5,x)-(Tk{p,pu,x)\~'- 1 dx) 2( ' 1 ^ds) 2 

Jo Jd 

(/ (/ \cr k {ps *Vs, {psus) *vs,x) - a k {p, {psus) *vs,x)\^dx) 2(1 ^ds )' 

Jo Jd 

(/ (/ Wk{p, {psus) * Vs,x) - a k {p,pu,x)\^dx) 2( ' 1 7) ds) 2 

J 0 J D 


< 


kfc|l- pim Er I {psus) *VS- P u \%(Ll) + I PS *VS- P\ 2 Ll(Ll) 


In the last line, we interpolated and used to the lipschitz Hypothesis 11.31 Let us now explain why these terms 
go to zero. Using Proposition 17.81 one can extract a subsequence such that p n —► p in L 2 ~ for all q < 2, P 
almost surely. Using the method on page 23 of [21] . one can use the renormalizations of the transport equation 
to upgrade this convergence to P almost surely in C t {Lf.~). Finally, appealing to Lemma [7/91 and Vitali we see 
that the quantity about tends to zero. 

This convergence suffices to complete our usual method and identify 


OO /» 

M t {4>) = V] / / pa k {p,pu,x)dxdfi k {s). 

k=1 J o Jd 
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We may also verify the energy bounds m using the lower-semicontinuity of the norm in the usual way. The 
desired continuity and measurability conditions imposed in Part [T| of Definition [L5] follows from the construction 
of the filtration L 0 - □ 

This completes the proof of our main result. 

Appendix A 

A.l Random Variables on Topological Spaces and the Skorohod Theorem 

Let (LI, F, P) be a probability space and (E , r, B T ) be a topological space endowed with its Borel sigma algebra. 
A mapping X : LI —> ( E , r) is called an ll E valued random variable” provided it is a measurable mapping 
between these spaces. Every E valued valued random variable induces a probability measure on ( E,t,B t ) by 
pushforward, which we denote Po A -1 . A sequence of probability measures {Pn}^! on B T is said to be “tight” 
provided that for each £ > 0 there exists a r compact set K^ such that F n (K^) > 1 — £ for all n > 1. 

A collection {X t }J =0 is an E valued stochastic process provided that for each t, X t is an E valued random 
variable. An E valued stochastic process is progressively measurable with respect to the filtration {.F t }^_ 0 
provided that for each t <T, 

X | [0 , t] :nx [0,t] -t ( E,t,B t ) 
is measurable with respect to the product sigma algebra F t x B([0,t]). 

Definition A.11. A topological space (E,t) is called a Jakubowski space provided there exists a countable 
sequence { G k }^ =1 : E —> R of t continuous functionals which separate points in E. 

Our main interest in such spaces is the following fundamental result: 

Theorem A.12. Let ( E,t ) be a Jakubowski space. Suppose that {Xk}k> l is a sequence of E valued random 

~ y. y. Y y. y. 'J OO 

variables on a sequence of probability spaces {(fA, F k , Pfc)}fc>i such that < P& o Xjf 1 l is tight. 

I J k—1 

Then there exists a new probability space P) endowed with an E valued random variable X and a 

sequence of “recovery” maps {Tk}ffLi 


T k :(fl,E,F)^(n k ,E k ,P k ) 


with the following two properties: 

1. For each k, the measure P^ may be recovered from P by pushing forward T k . 

2. The new sequence {Afc}fc>i := {X k oTiJim converges P a.s. to X(with respect to the topology t). 

Proof. This result is a combination of the versions of the Skorohod theorem proved in [16| and [26]. It can 
be proved by modifying the proof in [16] in a very slight way. Namely, at the point in the proof where the 
classical Skorohod theorem for metric spaces is applied, one may apply the Skorohod theorem in |26l to obtain 
the recovery maps. □ 

Remark A.13. It is straightforward to check that the following are examples of Jakubowski spaces: Polish 
spaces, dual spaces of separble Banach spaces B wt endowed with the weak star topology, and Ct(B w ) for reflexive 
Banach spaces B. 

Also, the class of Jakubowski spaces is closed under countable products. In particular, given a Jakubowski 
space ( E,t), E°° is also a Jakubowski space with respect to the r product topology. Similarly, for finite products 
of different Jakubowski spaces. 
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A.2 Series of One Dimensional Stochastic Integrals 

By a stochastic basis, we mean a probability space (fl, J 7 , P) together with a a filtration { and a collection 

{Pk}kL i °f one dimensional Brownian motions. 

Proposition A.14. Let (Q,J r ,W‘,{J 7t }J_ 0 ,{l3ky£L 1 ) be a stochastic basis endowed with a collection of {j rt }^ =0 
progressively measurable proceeses {f k }kLi '■ H x [0, T] — > K, such that 



E p [f k (s)} ds < 


oo. 


Then we may construct an {J rt }^l 0 martingale {Mt}J=Q with P a.s. continuous paths of the form 

oo „t 

M t = V / fk{s)dp k {s). 

k—1 J ° 

The series above converges uniformly in time in probability and the quadratic variation process is given by 


oo „t 

= fk{s,iw)ds 

k= i Jo 


Proof. This is a consequence of the Kolomogorov Three Series theorem and the construction of the one dimen¬ 
sional stochastic integral. See Krylov nn for more discussion. □ 

The next lemma, taken from [5], provides a procedure for identifying a continuous, adapted process as a 
series of one dimensional stochastic integrals. 

Lemma A.15. Let (11, T, P, {Pk}kLi) be a stochastic basis endowed with a continuous o mar¬ 
tingale Moreover, suppose the following are also {Ft}J = 0 martingales 

1. (uj,t) ^ fk(“, s ) ds 

2. (w, t) -»• Mj(w)/3 t fc (w) - / 0 ‘ fk(s)ds (for each k> l) 
then the process {Mt}J=$ may be identified as 

oo ,.t 

= V / fk{s)d(3k{s). 

k=i Jo 

A.3 The Space of Weakly Continuous Functions in Lff 

This section contains a useful tightness criterion for probability measures over the topological space C t {[LP] w ). 

Lemma A.16. Let {f n }((L i a sequence in Lf°(L™) with 1 < m < oo. Suppose that the following two 

criterion are met: 

1. sup„ < oo 

2. For all f> £ C(?°(D) in a dense subset of Lff, the following sequence in Ct is equicontinuous 

(*->• f fn{x,t)(j)(x)dx\ 

l JD J n=1 

Then there exists an f £ Ct{[Lff] w ) and a subsequence such that 

fn k t / in C t {[L™] w ). 
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Proof. See the Appendix in Lions [20j . □ 

A straightforward application of the lemma above yields 
Corollary A.17. For any positive M, integer k and q > 1, the following sets are compact in Ct([Ltf\ w ) 

{/ G C t ([LP] w ) | \f\ LT {L T) + \dtf\ L , (w -K n < M} 

The tightness criterion can now be stated as follows 

Lemma A.18. Let {/„}“Li be a collection of Ct([L™] w ) valued random variables, each defined on a probability 
space (p, n ,J r n ,F n ) such that 

1. 

SUp E P " | /n | L“ (L£*) < OO 

n 

2. For any <f> G C£°(D), there exists an integer k, 7 > 0, and p > f such that 

SUp E P " [\(fn(t) ~ fn(s),f))\ P ] < \(f\ P ck \t- S| 7P 
n x 

for all 0 < s,t < T. Then the sequence of induced measures {P n o are tight on Ct([L™] w ). 

Proof. Enumerate a countable collection in Cf°(D) which is dense in L™ . The second hypothesis of 

the lemma implies that for all s < 7 and j > 1 

supE|(/„,</>,•)|*Up < \(fj\ p ck . 

n> 1 * 

Choosing a > 0 sufficiently small to apply the Sobolev embedding theorem gives 

SU P E|(/ n , 4>j)\ P C a < \<t>j\ P ck . 
n> 1 * 

Given a small number f > 0, define a set K^ by 

OO 

*£ = {/e L?{L™) I \f\ L ? { L T ) < MC 1 } n fl {/ G L?{L™) I \{f,<t>j)\cf < (2 J T 1 )'l^lc*} • 

i=i 

Lemma IA. 161 implies this set is sequentially compact in C't([L r f l ] w ). Sequential compactness and compactness 
are equivalent in Ct([L™]tu). Applying Chebyshev, then using the uniform bounds, we find that 

supPo f-\KD<^ 


Appendix B 

B.4 Weak Convergence Upgrades 

The following lemma is simple, but fundamental enough to state explicitly. 

Lemma B.l. Let E,F be Banach spaces and use E w to denote the space endowed with its weak topology. Let 
T : E —>■ F be a bounded linear operator. Suppose the sequence {/n}^L 1 converges to f in Ct(E w ). Then 
{Tfnf^Li converges to Tf in Lj(F w ) for all 1 < q < 00. 

Lemma B.2. Make the same assumptions as in Lemma \B. 1\ above. In addition, assume T is compact. Then 
{Tf n }™ =1 converges to Tf in L'l(F) for all 1 < q < 00. 
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Proof. Since bounded operators preserve weak convergence, for each t £ [0,T] we have Tf n (t) Tf(t ) weakly 
in F. If T is compact then the convergence is strong. Combining the uniform bounds in Ct(E w ) with the Vitali 
convergence theorem gives both claims. □ 

Lemma B.3. Let (Ll,F,g) be a finite measure space. Let {f n }^Li in L p (Ll,F, p) converge weakly to f € 
L p (LI, T, p). Moreover, assume there is a convex function p : R —> R such that {p(fn)}%L i converges weakly to 
ip(f) in L 1 (LI, T, p). Denote by C the subset o/R where p is strictly convex. 

Then there is a full p measure set LI' such that converges pointwise to f(ui) for all uj £ C fl LI'. 


B.5 Some Tools from the Deterministic Compressible Theory 

The following result is a consequence of the Div Curl lemma. Denote 7 Zij = dij A -1 , understood to be well 
defined on compactly supported distributions. 

Lemma B.4. Let D be a smooth, bounded domain and rj a smooth cutoff. Let B be a Banach space. Suppose 
{fn}n Li converges to f in Ct([L p ] w ) and {g n }^L i converges to g in Ct([Lf] w ). Also, assume the embedding 
L r x ^f B is compact, where 2 + | = - < 1. 

Then the following convergence holds: 


V (fn^ij[V9n] ~ 9nKi,j[vfn]) T) (fKij[V9\ ~ 9^i,j[vf}) 
weakly in L r fi(B) for all 1 < m < oo. 

Proof. Combine the corresponding result in Feiresil 0 with Lemma |B.2f using the compact injection operator 
from L r x to B). □ 

Next we collect some properties of the Bogovoski operator B. Recall that classically, B\g\ is defined to be 
the solution to the problem 

fdiv„ = 9 in D (B 

I v = 0 on dD 

for g £ L P (D) such that f D gdx = 0. For our purposes, it is useful to have an extension of this operator to the 
negative Sobolev spaces. We recall a result from [15]. Define the space Lq(D) = {/ £ L P (D) \ f D fdx = 0}. 

For s £ [0,1] define W S ’ P (D) := W S ’ P (D) D Lg(D). Furthermore, let W~ S ’ P (D) := [W S ’ P '(D)}'. 

Theorem B.5. Let 1 < p < oo and s £ [—1,1]. Then there exists a bounded linear operator B : W S,P (D ) —> 
\W s+1 ’ p (D)] d such that div£%] = g for all g £ W S ’ P (D). 


B.6 Lemmas on Parabolic Equations 


The following lemma provides an energy equality for sufficiently integrable weak solutions to the parabolic 
Neumann problem driven by a rough velocity field. 


Lemma B.6. Let u £ L^fWf’ 2 ) and p > d. Suppose p £ Lf°(L p ) is a distributional solution of 


dtp — eA p + di v(pu) = 0 


fip- = 0 

dn U 


.p(z,0) = po{x) 

Then for all times 0 < t < T, the energy identity holds: 


in D x [0, T) 
in dD x [0, T] 
in D 


(B.2) 


f p 2 (t)dx + e f f \X/p\ 2 dxdt = 7 - f p^dx — ^ / f div u p 2 dxdt. 
' D JO J D 2 J D 2 J Q J D 


We also need a variant of the usual Ll(Wf ,p ) estimates for the parabolic Neumann problem. A similar result 
is proved in the appendix to [3j. The lemma below states that by giving up the optimal exponent, one can 
retain a form of the usual estimate even if the solution. Recalling the splitting from section [3] Choose r such 
that 2] is an even integer and let tk = kr for k = 0,..., -y. Define h^ et by (13.11) . 
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Lemma B.7. Let 1 < p < oo and F € Lf x . Suppose that p solves 


Then for all q < p 


dtp - eAp = F 

in 

d t p = 0 

in 

<>£. — 0 
dn U 

on 

p(0,x) = p 0 

in 


D X UfcLo {t2k,t2k+l] 

D x UfcL 0 (t2k+l,t2k+2] 

dD x [0, T] 

D 


+ \P\L q t (Wf’ q ) 


< C(e, d) (\Po\ w ^p + ■ 


(B.3) 


(B.4) 
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